arXiv: 1501.00215v3 [quant-ph] 24 Aug 2015 


Few-Body Systems manuscript No. 

(will be inserted by the editor) 


N.L. Harshman 


One-Dimensional Traps, Two-Body Interactions, 
Few-Body Symmetries 


I. One, Two, and Three Particles 


Received: date / Accepted: date 


Abstract This is the first in a pair of articles that classify the configuration space and kinematic 
symmetry groups for N identical particles in one-dimensional traps experiencing Galilean-invariant 
two-body interactions. These symmetries explain degeneracies in the few-body spectrum and demon¬ 
strate how tuning the trap shape and the particle interactions can manipulate these degeneracies. 
The additional symmetries that emerge in the non-interacting limit and in the unitary limit of an in¬ 
finitely strong contact interaction are sufficient to algebraically solve for the spectrum and degeneracy 
in terms of the one-particle observables. Symmetry also determines the degree to which the algebraic 
expressions for energy level shifts by weak interactions or nearly-unitary interactions are universal, 
i.e. independent of trap shape and details of the interaction. Identical fermions and bosons with and 
without spin are considered. This article sequentially analyzes the symmetries of one, two and three 
particles in asymmetric, symmetric, and harmonic traps; the sequel article treats the N particle case. 

N.L. Harshman 

Department of Physics 

American University 

4400 Massachusetts Ave. NW 

Washington, DC 20016-8058 USA 

Tel.: -f 1-202-885-3479 

Fax: -t 1-202-885-2723 


E-mail: harshman@american.edu 




2 
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1 Introduction to Part I 

The focus of this pair of articles is the non-relativistic, one-dimensional, few-body Hamiltonian with 
the following characteristics: (1) Each particle has the same mass and experiences the same trapping 
potential. (2) There is a two-body interaction term for each pair that depends only on the distance 
between particles. (3) Each particle has a finite number of internal levels that do not participate directly 
in the trap or two-body interactions. The particles could be distinguishable, or they could be identical 
bosons or fermions. The total Hamiltonian for the system can be expressed as 

N N 

(la) 

i=l i<j 

Denoting each canonical pair of particle observables by [Qj, Pk] = iSjk and choosing natural units, the 
one-body Hamiltonian for particle i is 

Hi = IP.^ + V\Q,). (lb) 

The two-body interaction term Vij has the Galilean invariance property Vij = V^dQi — Qjl). Particular 
attention is focused on the contact interaction, expressed in particle coordinates q = (gi, <72,. •., Qn) 
as 

qjl) = 9j)- (Ic) 

The goal of this pair of articles is to classify the symmetries of the few-body Hamiltonian for 
the cases of no interaction, general interaction, and unitary limit of contact interaction and then to 
demonstrate how these symmetries can be used to calculate spectral properties and understand univer¬ 
sal features. Two classes of symmetries are considered: configuration space symmetries and kinematic 
symmetries. By configuration space symmetry, I mean the group of transformations of conhguration 
space ~ that are represented as unitary operators that commute with . Configuration 
space symmetry includes the permutation group of identical particles, but it also can include parity or 
emergent symmetries depending on the trap and interaction potentials. Kinematic symmetry 
is realized by the group of all unitary operators that commute with . The kinematic symmetry 
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group necessarily contains the configuration space symmetry and time translation as subgroups. A key 
insight is that the dimensions of the irreducible representations of the kinematic symmetry group (if 
properly identified) explain the degeneracies in the spectrum of . 

This first article analyzes the symmetries of one, two, and three particles. The configuration space 
symmetries and kinematics symmetries are developed incrementally, and the ways in which the trap 
shape and the two-body interaction effect the symmetry are explained with examples. For systems 
with few degrees of freedom, the order of finite symmetry groups are small, so explicit calculations and 
applications are included. Additionally, the symmetries of one, two and three particles can be visualized 
using familiar geometrical methods and analogies. The sequel article treats the general case of N 
particles. In that case, the formal, algebraic machinery of group representation theory demonstrates 
its power. However, the price is a higher degree of abstraction and the necessity of computer-based 
algebraic methods. 


1.1 Motivation 

The model Hamiltonian Q has a long history inspired by applications to atomic, molecular, nuclear and 
condensed matter physics. Going back to the beginnings of quantum mechanics, various subfields have 
given different names (e.g. Stoner Hamiltonian, Tonks-Girardeau gas, Lieb-Liniger model, no-core shell 
model) to particular instances of the model and its higher dimensional generalizations. There is a large 
mathematical physics literature on the one-dimensional model, and certain cases of are exemplars 
of solvability in few-body and many-body systems [mil nig. The increasingly precise preparation, 
control and measurement of ultracold trapped atomic systems in effectively one-dimensional traps |6] 
is driving another surge of theoretical interest in this few-body model, c. f. Illllllllllllllllll 

lIllliailSlinilllllHlEolEIlElEllMlEaEaEaEQlEgEgEIllslEllMllaaiMlEaislIlo]. 

Few body properties can drive the dynamics of many-body cold atom trapped systems, like trap loss 
and equilibration, and few-body observables may be more directly accessed by tunneling rates and 


spectroscopic methods. 
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Although group theory has a long history of being productive in quantum mechanics, the “Grup- 
penpest’l^ can be so frustrating that it is customary to begin with an explanation of why all this 
mathematical apparatus is worth the effort. The essential claim is that the symmetry classifications 
provided in this article can be exploited for qualitative, analytic and numerical studies of few-body sys¬ 
tems trapped in one dimension and they provide a unifying framework for this recent wave of analysis. 
These methods solve or simplify numerous questions about the spectrum, degeneracy and dynamics, 
including the following: 

— identical particle symmetrization, 

— perturbation theory from the non-interacting to the weak interaction limit, 

— perturbation theory from the unitary limit of the contact interaction to the nearly-unitary limit, 

— methods of exact diagonalization in truncated Hilbert spaces, 

— perturbation theory for not-quite identical particles, 

— adiabatic or non-adiabatic particle dynamics under variation of interaction parameters or trap 
shape, and 

— trial wave functions for variational or Monte Carlo methods. 

As a preview of the kind of results that symmetry classification and calculations provide, see Fig. Bit 
depicts how level splitting in the near-unitary limit of the contact interaction depends on trap shape 
for four particles. Depending on whether the particles are fermions or bosons, with or without spin, 
only certain energy levels can be populated. The method of calculation is developed later in the paper, 
but the main idea is that near unitarity, level splitting is determined by the tunneling amplitudes 
of adjacent particles and these tunneling amplitudes depend on the shape of the trap. The energy 
eigenstates can be found by diagonalizing a tunneling operator, and these eigenstates carry irreducible 
representations of the symmetric group for four particles S 4 and for the parity symmetry. 

Many applications of the representation theory of the symmetric group already exist in the recent 
literature; a few examples relevant to these articles are (nmaiiiiiiiiEsiiMiiM!. Parity is also widely 
exploited, and the special symmetries of harmonic traps are often explicitly or implicitly invoked. 
The focus of this article is to see how much more solvability is provided by additional configuration 
^ See the Introduction to m for a discussion of the Gruppenpest. 
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Fig. 1 Level splitting diagram for four particles in three symmetric traps with contact interactions: (a) double 
well; (b) infinite square well; (c) V-shaped or cusped well, i.e. softer than harmonic. The thick band on the 
right in each figure is the 24-fold degenerate ground state energy level for four distinguishable particles in the 
unitary limit (/ —>■ oo of the contact interaction. This level carries two representations of S 4 simultaneously; one 
copy corresponds to particle permutation symmetry P4 and the other to ordering permutation symmetry. In 
the near unitary limit, ordering permutation symmetry D 4 is broken by tunneling and the energy levels split 
into irreducible representations of P4 with either even or odd parity. The thickness of the line indicates the 
degeneracy of these levels for distinguishable identical particles without spin. For example, the irrep labeled [4]^ 
is the non-degenerate, positive parity, totally symmetric spatial state. It can be occupied by bosons with any 
number of internal levels, or fermions with at least four internal levels. The three trap shapes are distinguished 
by the ratio of tunneling amplitudes tju, where t is the tunneling amplitude for the left-most or right-most 
particle to exchange with the adjacent inner particle and u is the tunneling amplitude for the two inner particles 
to exchange. To first order, both are proportional to l/g. The following ratios have been chosen to illustrate 
the trap dependence of these amplitudes: (a) t/u = 2.9; (b) t/u = 1; (c) t/u = 0.3. The idea is that (a) for 
double wells tunneling in the middle is suppressed so t > u; (b) for inhnite square well the potential is uniform 
inside the trap and so (for low particle density) t and u are approximately the same; (c) for softer wells, there 
is more phase space in the middle of the well so u > t. For harmonic wells, t/u « 0.762 (c.f. |29l 1301 [371139p . 
In subfigure (c) a more extreme ratio is depicted, corresponding to a V-shaped or cusped trap. 


space and kinematic symmetries inherited from the trap shape and the Galilean invariance of the 
interactions. We know that in the case of the infinite square well and contact interactions of any 
strength, there is enough symmetry to provide integrability, i.e. the Bethe ansatz solutions (c.f. [Hi). 
The experimental tunability of few-body symmetries and the close connection between finite groups 
and integrability |43| |44j suggest novel possibilities for embodying mathematical structures in ultracold 
atomic systems. 










































6 


Symmetry also aids the study of “universal” few body phenomena, a term used (with some local 
variation) to describe dynamical effects that do not depend strongly on the particular details of the 
constituent few body systems or on the nature of their interactions. See Hisi sg for discussions 
of universality in one-dimension. Universal properties established in atomic systems can also reveal 
themselves in few-body systems at the chemical or nuclear scale. Universality can drive the dynamics 
of coherence, entanglement and equilibration in certain many-body systems. One approach to univer¬ 
sality is to figure out how much about the few-body system can be inferred from the symmetries of 
without specihc knowledge of the trap or the interaction. The relationships among trap shape, 
interaction, and permutation symmetry determine which properties of the system can be algebraically 
solved for in terms of one-particle observables. The degree to which a few-body system possesses this 
kind of ‘algebraic solvability’ is at least some component of universality. The best example is provided 
by the unitary limit of the contact interaction, which has enough symmetry to be exactly solved for 
any N given the one-particle spectrum [51 UHl [HI HSl IHl ttH [13 [13 [HI EH El] . The question of how 
level splitting in the weak interaction limit and near-unitary limit depends on trap shape is a theme 
that runs throughout this pair of articles. 

Symmetry methods also provide geometrical insight into the highly-abstract interplay of trap shape, 
interaction, spin, and particle symmetrization. Especially for low particle numbers, symmetries can be 
pictured and manipulated in the mind. To a large extent, the geometrical constructions and geometrical 
methods applied in the works m uni E3 EOl EH El 113 motivated this article. In this first article, I 
argue that by analyzing the cases of two and three particles geometrically, we get insights that can 
guide us for higher particle numbers where more abstract methods are required. 


1.2 Outline of the Articles 

The next section of this article explains which configuration space and kinematic symmetries are pos¬ 
sible for one particle in asymmetric and symmetric traps, and explains the extra kinematic symmetry 
that occurs for the harmonic trap. The third and fourth sections develop symmetry classifications and 
techniques for two and three particles. In each scenario, the non-interacting case is considered first, 
then the interacting case (including weak interactions), and finally the unitary limit of the contact 
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interaction (including the near-unitary limit). For three particles, state permutation symmetry and 
ordering permutation symmetry are introduced as useful concepts complementary to the more familiar 
particle permutation symmetry. Along the way, a variety of applications, diagrams and figures are 
included that attempt to make the symmetry methods more concrete and less abstract. This article 
ends with a conclusion that reflects on what this symmetry analysis says about universality in this 
model. 

The second article in this series derives the general form of the symmetry classifications for N 
particles. It is necessarily more technical (and has fewer pictures). After an introduction that gives 
the expressions for the minimal configuration space and kinematic symmetries inherited by the con¬ 
struction of the few-body system from the one-body systems with two-body interactions, the next 
section gives an overview of the symmetric group Sn and its representations. Definitions, notation and 
conventions necessary to extend these methods for N > 3 are briefly reviewed. In particular, a kind of 
Sat representation space called a permutation module is shown to be especially useful for the analysis 
of N identical particles. The third section establishes the symmetries for N non-interacting particles 
and describes the geometric realization of particle permutations and other symmetries in configura¬ 
tion space. The irreducible representations for the minimal kinematic symmetry group are derived 
and state permutation symmetry is used to construct a complete set of commuting observables that 
facilitates identical particle symmetrization. A final result of this section establishes the isomorphism 
between the bosonic non-interacting spectrum and the fermionic spectrum (which remains invariant 
under contact interactions). The fourth section classifies the symmetries for N particles interacting 
via two-body Galilean invariant potentials. The symmetries of two-body matrix elements are derived 
and state permutation symmetry makes another appearance, this time as a property of the contact 
interaction. The two body matrix elements are used to analyze level splitting in the weak interaction 
limit and I conjecture that algebraic solvability is lost for more than five multicomponent particles. 
The fifth section treats the unitary limit of the contact interaction. Ordering permutation symmetry 
emerges as new symmetry of the system, and the near-unitary limit can be understood in terms of 
symmetry breaking of ordering permutation symmetry by tunneling among different sectors of con¬ 
figuration space. The final and concluding section of both articles summarizes how the main results 



relate to the question of universality and describes some possible further extensions and applications 
of this work. 

1.3 A Few Notes about Group Notation 

These articles are addressed to several distinct audiences, including novices and experts, interested 
in low-dimensional, trapped ultracold atomic systems, general quantum few-body systems, and/or 
mathematical physics. I have attempted to clearly signpost the content into sections, subsections and 
subsubsections so that readers can pick and choose what matches their interests and background. The 
first article is more pedagogical and less technical. The second article presumes more familiarity with 
group representation theory, but most necessary ideas are developed in this hrst article. 

Another challenge when providing clarity for a diverse audience is in the choice of notation. This 
is particularly important when taking about symmetry, groups, and group representations because 
different physical symmetries may be isomorphic to the same abstract group, and groups usually have 
multiple inequivalent representations. The next few subsubsections provide a brief introduction to the 
notations for symmetries, abstract groups, and their representations that will be used in these articles. 

1.3.1 Groups and Representations 

The configuration space symmetry group for an N particle system with Hamiltonian Q is denoted Cat. 
For the particular case when there are no two-particle interactions, the configuration space symmetry 
is denoted and for the unitary limit of the contact interaction the group is denoted CJ/. The 
kinematic symmetries are similarly denoted K^, and K“. 

These symmetry groups are isomorphic to abstract groups. The specific abstract group depends on 
the shape of the trap. For example, for non-interacting particles in an harmonic trap, is isomorphic 
to U(iV), the abstract group realized by unitary N x N matrices. To highlight the distinction between 
isomorphic and equality, I write ^ U(A^). 

For a given group G, up to three different representations are in the analysis of this article: 

1. Unitary irreducible representations, or irreps. Almost all groups in this article are hnite or compact, 
and they have a finite or countable number of finite-dimensional irreps. Other representations are 
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built out of direct sums of irreps. The labels or notations for irreps depend on the group. As an 
example, pretend the symbol 4|k labels a particular irrep of the group G. The dimension of the irrep 
is d{G; 4|l), or (i(4|l) if the group is obvious from the irrep label. The d(4|l)-dimensional unitary matrix 
representation of 5 G G is denoted (g). The complex vector space that carries the representation 
is ~ 

2. Hilbert space representation. Every element of a symmetry group g € G is represented by a unitary 

operator on the Hilbert space, denoted U(g) or g. The Hilbert space can be decomposed into irreps 
of the symmetry group. For example, if the irreps of G are 4|k, <C>, and Jit, then 

H = H* 

Each subspace of H could be a single irrep, i.e. ^ , but generally the subspace is a 

tower of irrep spaces 

i 

where i is a label or set of labels that distinguish different copies Hf ~ of equivalent irreps of 
G in 

3. Configuration space representation. This refers to the action of Cat (or the special cases C% or C“) 
on the configuration space {qi,q 2 , ■ • ■, (Zat) S = K'^. The representation of g S Cat is denoted 
0{g). Typically, this representation of Cat is not irreducible. 

1.3.2 Symmetric Group 

One obvious symmetry of the model Hamiltonian Q is the group of particle permutations Pat. The 
configuration space symmetries Cat, C%, and all must contain Pat as a subgroup. The elements 
p € Pa/ can be described either in permutation notation or cycle notation. For example, the same 
permutation p can be written either as permutation {312} or three-cycle (132). Both notations for p 
describe the map in which particle 1 is replaced by particle 3, particle 2 is replaced by 1, and particle 
3 is replaced by 2. 

The group Pat is isomorphic to the abstract group Sat, the symmetric group on N objects. Two 
other groups described in later sections are also isomorphic to symmetric groups, the group of state 
permutations ‘^yu\ on the state composition lv\ and the group of ordering permutations On on N 
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particles. The symmetric group S^v has order N\ and it has irreps labeled by whole number partitions 
of N denoted [fj] = [ 121^2 ■ • • /^r] where These irreps are sometimes depicted by Ferrers 

diagrams (also called Young diagrams), which are r rows of boxes with /i^ boxes in each row. The Sjv 
irrep spaces are denoted Ai and the matrix representation of p S Sat on is denoted (p). 

A few notes and examples with Sjy for iV = 1, 2, and 3: 

1. The group Si ~ Zi is trivial. 

2. The group S 2 Z 2 has two elements e and (12) and it is abelian. It has two one-dimensional irreps 

labeled by the Ferrers diagrams ^ and B, or more compactly by partitions [2] and [11] = [1^]. 
For the trivial, symmetric representation, we have = 1 and for the faithful, antisymmetric 

representation ^ 12 ) = — 1 . 

3. The group S 3 has six elements in three classes: the identity e, three two-cycles (or transpositions) 
(12), (23), and (31), and two three-cycles (123) and (132). This group is not abelian, and since 
it is not abelian, no faithful representations can be one-dimensional. There are three irreps: [3] or 
rm (one-dimensional, symmetric); [ 21 ] or P (two-dimensional, faithful); and [ 111 ] = [ 1 ^] or ^ 
(one-dimensional, antisymmetric). 


1.3.3 Point Groups 

Most of the configuration space symmetries considered in this article are point groups. Point groups 
are orthogonal transformations of the iV-particle configuration space = M.^, and the set of all 
possible point groups for a given dimension is completely characterized [48] . When many people think 
of symmetry, it is the geometrical realization of point-group invariant objects that they envision. One 
important class of point groups are the finite Coxeter groups. These are generated by reflections in 
Euclidean space, they are the symmetries of regular polyhedra, and their categorization is closely 
related to the structure of simple Lie algebras. 

The maximal point group for is the group 0{N) of all orthogonal transformations in N dimen¬ 
sions, i.e. all reflections and rotations. All other point groups in N dimensions are subgroups of 0{N). 
Here are a few facts about point groups in low dimensions useful for understanding this article: 
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1. In one dimension, there are only two point groups. One is the trivial group that contains just the 
identity e. The other is the group 0(1) that contains the identity e and a single reflection 77. These 
groups are isomorphic to the abstract cyclic groups Zi and Z 2 , respectively. Since these groups are 
both abelian, all irreps are one dimensional. Irreps of 0(1) are labeled by tt = ±1. 

2. In two dimensions, besides 0(2), there are two series of finite-order point groups. This article em¬ 
ploys several of the dihedral groups Dj,, finite groups of order 2k that include k rotations (including 
the identity) and k reflections. The group Di ~ Z 2 is the symmetry of a butterfly, the group D 2 is 
the symmetry of a rectangle, the group D 4 is the symmetry of a square. For k > 2, the group 

is not abelian and so its faithful irreps are not one-dimensional. 

3. In three dimensions, besides 0(3), there are seven infinite series of finite-order point groups, seven 
other finite-order point groups, and four other continuous point groups. These groups are familiar 
to some from chemical or solid states physics; see [50] for a palatable introduction to these groups 
and their irreps. There are multiple conventions for the notation of three-dimensional groups and 
irreps (Schonflies, Coxeter, orbifold, etc.); specific notations are introduced as necessary. 


1.3.4 ^ Fs-w Other Groups and Notes 

The abstract group of translation by a single parameter a; S M is denoted T^,. Irreps of T^, are one¬ 
dimensional and abelian and characterized by a single real number. Specific examples include: the 
group of time translation Tj represented on the Hilbert space by U{t) = exp(—iHt) with irrep labels 
called energy; and the group of space translations T^ represented by U(a) = exp(—iPa) with irrep 
labels called momentum. 

Finally, note the following: 

- Groups are generally denoted by capital Roman letters, e.g. Sat, Cat, U(7V), etc. 

— Vector spaces are denoted by capital calligraphic letters, e.g. the total Hilbert space 77, the spa¬ 
tial Hilbert space /C, the spin (or internal component) Hilbert space S, irrep spaces N4, or the 
configuration space Q. 
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— Operators on the infinite-dimensional Hilbert space have ‘hats’ like Qi and Ri. Their eigenvalues 
are usually lowercase like qi and Matrix operators on finite-dimensional spaces like irrep spaces 
M or configuration space Q are underlined, e.g. the representations D and O. 

— Ordered sequences of numbers or symbols are denoted by angle brackets, e.g. {v) = (3,0, 2,4,4,0) 
and (/i) = {a,P,a). Compositions of unordered numbers or symbols are denoted by floor brackets. 
For example, the compositions of the previous two sequences are = [002344J = [0^234^J 
and [/rj = [aa/?] = [a^/3J. The set of all sequences with a composition {vl is the same as all 
permutations of the composition, denoted P . The shape of a composition is the pattern of 
degeneracies in a composition, e.g. [v] = [2211] = [2^1^] and [/x] = [21], and always corresponds to 
a partition of the length of the sequence. 

— The non-negative integers {0,1, 2, 3,4, • • • } are denoted N. 

2 One-Particle Symmetries 

This section describes the configuration space symmetry group Ci and the kinematic symmetry group 
Ki for one-particle in asymmetric, symmetric and harmonic traps. The symmetries Ci and Ki are 
built from basic abstract groups that have only one-dimensional representations. These one-particle 
symmetry groups are the building blocks of the multi-particle analysis. 

Consider one particle in a one-dimensional trap and denote its spatial Hilbert space K,. The total 
Hilbert space % = JC®S is the tensor product of the spatial Hilbert space and the spin Hilbert space 
(discussed at the end of this section). All one-dimensional systems have at least the symmetry group 
T( of time translations. Although this observation seems trivial, this symmetry is enough to guarantee 
integrability for any one-dimensional system. The abelian, one-parameter group of time translations 
T( has one-dimensional irreps labeled by the energy e and the set of allowed energies determined by 
the Hamiltonian is the spectrum cti. Time translation group is represented by exponentiation of 
the Hamiltonian tj{t) = exp(— 

For a single particle trapped in one dimension, the energy spectrum ai = {co; ei) • • ■} is discrete, 
count ably-infinite and non-degenerate. An energy spectrum with this simple form excludes the impor¬ 
tant idealized case of infinite lattices and periodic boundary conditions. Further, a discrete spectrum 
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is only a low-energy approximation for wells with finite depth because it does not have a continuous 
piece. There are probably other interesting pathological cases not covered, however this kind of spec¬ 
trum does include double-wells, multiple-wells and all the greatest hits of one-dimensional solvability 
like the harmonic well, infinite square well, Polsch-Teller potential, Morse potential, etc. 

Eigenstates of are denoted by kets containing the spectral index 

H^\n)=en\n). (2) 

and the corresponding wave functions are 

(l>niq) = {q\n). (3) 

No functional dependence of e„ on n is implied, although algebraic or transcendental expressions 
certainly exist for specific solvable potentials. For convenience, sometimes the one-particle eigenstates 
will be denoted by state labels |a), |/3), jy), etc. with wave functions 4>aiq) and (for symmetric traps) 
parities tTq. 


2.1 Configuration Space Symmetries for One Particle 

The configuration space symmetry Ci is the group of all transformations of = M realized by 
operators that commute with the one-particle Hamiltonian H^. For an asymmetric trap, no such 
operators exist and Ci ~ Zi is the trivial group. 

For a symmetric trap, there is a single point about which reflections are a symmetry and Ci ~ 0(1) 
is the parity group. For symmetric one-dimensional wells, the quantum number n also determines the 
parity 

n\n) = {-ir\n). (4) 

Although not a trap (and outside the purview of this article), for a constant potential (e.g. no 
potential) the group Ci is the Euclidean group in one dimension Ei ^ 0(1) x T^, where T^ ~ is 
the group of spatial translations in and x denotes the semidirect product. See [49] for a discussion 
of symmetries and partial symmetries of lattice-like multi-well potentials. 
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2.2 Kinematic Symmetries for One Particle 

The one-particle kinematic symmetry group Ki is the group of all unitary symmetry operators that 
commute with H^, and therefore necessarily contains Ci. 

Asymmetric Traps: For asymmetric traps, the only symmetry is time translation so Ki ^ Tj. The 
irreps are one-dimensional, consistent with the non-degeneracy of ui and are labeled by the energy e„ 
or quantum number n. The spatial Hilbert space can be decomposed into irreps of Ki: 

OO 

/C = 0/C". (5) 

n—0 

Each summand KT is the one-dimensional irrep of time translation where time evolution is represented 
as iy^{t) = exp(— 

Symmetric Traps: For symmetric traps, the kinematic group is Ki ~ 0(1) x T^. Irreps are still 
one-dimensional and labeled by n. The decomposition of JC into Ki irreps is the same as ([^, except 
now parity 7r„ = (—I)" is also a good quantum number. Therefore the spatial Hilbert space K. also has 
a decomposition into sectors of fixed parity JC = /C’*' © IC~ where 

OO OO 

/C+ = and JC- = 0/C2'=+b 

fc =0 fc =0 

Harmonic Traps: For harmonic traps Ki ^ U(l) x Tj. Here U(l) is the group of transformations 
that changes the phase of the ladder operators a and d^. Define a unitary representation of U(l) by 
operators for (j) G [0, 27r) such that 

b = U{(j))aU'^ {(p) = exp(*(/))d. 

This transformation leaves = a^a +1/2 invariant and can be thought of as rotations in two- 
dimensional phase space. Again, there is no change to the decomposition of 1C into Ki irreps ([^. 

Note that for all three kind of traps, Ci and Ki are abelian groups. For abelian groups, irreducible 
representations are one-dimensional, and this is consistent with the assumption of a non-degenerate, 
discrete one-particle spectrum cri. 
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2.3 Including Spin 

Finally, if the single-particle Hamiltonian is spin independent and there are J spin components, 
then there is also a factor group of U(J) to the kinematic symmetry. The total Hilbert space for the 
one particle system is the tensor product of the spatial Hilbert space and the spin Hilbert space. 

-H =/C(g)5-L2(M)(g)C'^. (6) 

Any unitary operator that acts only on 5 ~ C'^ certainly commutes with . Further, if the internal 
components really are spin components of a particle with spin s, then J = 2s-\-l and the spin operators 
and Sz form a complete set of commuting operators for S that commute with the Hamiltonian. 


3 Two-Particle Systems 

The purpose of this section is to classify the types of symmetries found for two trapped particles in 
the case of no interaction, a general two-body interaction, and the contact interaction. In some sense, 
symmetry analysis does not provide anything remarkable or new for two particles. However, it provides 
a training ground for intuition about symmetries in a familiar setting and it is useful for contrast with 
more complex scenarios. Also, techniques and notation are introduced here that are be extended to 
the three particle case in the next section, and then to the fV-particle case in the sequel article. 

One of the most important ideas of this section that the degeneracy of the two-particle spectrum 
fT 2 can be explained by looking at the kinematic symmetry group K 2 . The dimensions of K 2 -irreps 
should be the same as the degeneracies in 172 ■ If not, that could signal the presence of an emergent 
two-particle symmetry, i.e. a symmetry that cannot be generated from one-particle symmetries and 
particle permutations. 


3.1 Two Non-Interacting Particles 

Consider the total non-interacting Hamiltonian constructed from the sum of one-particle Hamiltonians 


Hi = Hi + Hi 


i^[Pi+P2)+V\Q^) + V\Q^). 


(7) 
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Fig. 2 On the left, each pair of one-particle quantum numbers numbers is a composition \y\ corresponding to 
an energy in the two-particle non-interacting spectrum a®. The superscript denotes that the compositions 
like \ a^\ with shape [u] = [2] are non-degenerate and compositions like [a/3} with shape [v\ = [1^] are two-fold 
degenerate. There is a totally symmetric spatial state in both shapes of compositions, but only compositions like 
[af3\ contain totally antisymmetric states. Note that there is a one-to-one map (depicted with dotted arrows) 
from all compositions levels in cr® (on left) to all mixed compositions (a/lj (on right), and the map preserves 
the partial ordering. This is one way to depict the famous boson-fermion mapping for two strongly-interacting 
particles in one dimension. 

The two-particle, non-interacting spectrum, denoted ct®, remains discrete and countably-infinite, but 
unlike the one particle spectrum ai it is necessarily degenerate. Every energy € <^2 is associated 

to (at least) one composition [a(3\ of two energies G cri. Unless the specific values of one- 

particle energies are known, only a partial ordering of tr® i® possible. The lowest two energies in (T 2 
are unambiguous: E^qq^ = 2eo and = eg -|- ei. However, the comparison of Eyny and EyQ 2 y is not 

possible without specific knowledge of the values for eg, ei and 62 - For example, consider the potential 
= \<lY■ For z = 2 the spectrum is harmonic with Eyny = if[g 2 j, for 0 < z < 2 the spectrum is 
softer than harmonic with Eyny > if|^g2j, and for z > 2 the spectrum is harder than harmonic with 
U^iij < EyQ 2 \ - See Fig. [^for a depiction of the partial ordering that can be put on cr^ without specific 
knowledge of cti . 



3.1.1 Two Non-Interacting Particles: Configuration Space Symmetries 

The configuration space for a system of two particles is Fig. [^depicts equipotentials for 

six traps, and without interactions this two-particle system is equivalent to one particle navigating 
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Fig. 3 These are the equipotentials for two particles in a one-dimensional trap with (a) V{q) — (cusped); 
(b) V{q) — |g| (V-shaped); (c) V{q) — q^ (harmonic); (d) V{q) = (approximately hard wall); (e) V{q) = 
q* — 2q^ (symmetric double well); and (f) V{q) = q'^ — IfSq^ — 2q^ (asymmetric double well). The horizontal axis 
is the qi axis and the vertical axis is q 2 ■ These figures all have particle exchange symmetry, which is realized in 
configuration space as a reflection across the line qi = q 2 - Antisymmetric spatial states must have a node on this 
line. Potentials (a)-(e) are also symmetric under spatial inversion of each particle individually, corresponding 
to horizontal and vertical reflections. For (a), (b), (d) and (e), the total symmetry group, combining parity and 
particle exchange is isomorphic to the non-abelian, two-dimensional point group D4. This order 8 group is the 
symmetries of a square and includes the identity, rotation by ±7r/2, rotation by tt, and reflection across three 
axes. The equipotentials of the harmonic well have the maximal point symmetry in two-dimensions: 0(2), i.e. 
rotations by any angle and reflections across any axis. 

these two-dimensional potentials. At a minimum, the configuration space symmetry group C® for two 
identical, non-interacting particles must contain as subgroups two copies of the one-particle symmetry 
group Cl. The particle permutation group P 2 ~ S 2 must also be a subgroup. This subsubsection 
establishes that the right way to combine these symmetries is 

D P2 >< Cf 2. ( 8 ) 

Asymmetric Trap: The case with the absolutely minimum symmetry possible is the asymmetric 
trap. Then Ci ~ Zi. and Cq ~ S 2 . Particle exchange (12) acts on by 


Q{12){qi,q2) = ( 92 , gi). 
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This representation on is isomorphic two-dimensional point group denoted Di, the dihedral group 
generated by single reflection along the line qi = 92 - 

The spatial Hilbert space can be decomposed into subspaces corresponding to the irreps of S 2 : 

= (9) 

where each of is a tower of symmetric states and is a tower of antisymmetric states. The 
Hilbert space representation of particle exchange is the unitary operator 17(12): 

17(12)|a/3) = |/3a). (10) 

Symmetric basis vectors \aa) = | m ) are invariant under 17(12) and are elements of the irrep tower 
However, the particle basis energy eigenstates |q;/ 3) and \j3a) {a ^ /3) do not belong to irrep 
towers. Instead, define the following simultaneous eigenvectors of Hq and 17(12): 

Iffl) = (|a/3) + \l3a)) G 

i>^^(|a/3)-|/?a))G/Cli\ (11) 

Symmetric Trap: For a spatially symmetric trap, the one-particle configuration space group is parity 
Cl ~ 0(1). Without interactions, each particle can be independently spatially inverted. Denote each 
particle’s inversion operator on by Ui such that 711(51,(72) = (—51,52) and 712(51,52) = (51,-52). 
These symmetries are represented in the Hilbert space on the particle basis as 

fliWi n 2 ) = (-l)”'|ni 712 ). ( 12 ) 

Including these two operations, the symmetry group C 2 for a symmetric trap has eight elements: 

e, (12), 77i, 772 , (12)77i, ( 12 ) 772 ,7Ti 772 , and ( 12 ) 77 i 772 . (13) 

This group is isomorphic to the point group of a square D 4 . The corresponding transformations on 
are 


77(0), 77 ( 71 / 4 ), 77(7r/2), 77(0), 77(7r/2), 7?(-7r/2), R{-k), and 77(-7r/4), 


(14) 
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where -R( 0 ) is a rotation about the origin by (j) and is a reflection across the line making an 

angle (f) with the qi axis. See Fig. and contrast the first five subfigures, which all have at least D 4 
symmetrjQ while the last subfigure only has Di symmetry. 

Note that the group D 4 is not abelian, e.g. (12)iTi = 772 ( 12 ). Therefore D 4 is not isomorphic to 
the direct product S 2 x 0(1) x 0(1) which would be abelian. Instead it is isomorphic to 

D 4 - S 2 X (0(1) X 0(1)) = S 2 X 0(1)^2_ (15) 


The notation x stands for the semi-direct product and captures the fact that the particle exchange 
(12) S S 2 conjugates elements 77i and 7 T 2 and therefore acts as an automorphism of the normal, abelian 
subgroup 0(1)^^. 

The group D 4 has five irreducible representations ED], four unfaithful one-dimensional irreps de¬ 
noted Ai, A 2 , 7?i and B 2 and and one faithful two-dimensional irrep denoted E. The spatial Hilbert 
space can therefore be decomposed like 

/C = ©/C^i ©/C®^ ©/C^. (16) 


The first four irrep towers in (16 1 contain one-dimensional irreps that have positive total parity 77 = 
77 i 772 and the two-dimensional irrep E has negative total parity. As an example, the energy levels 
included in Fig. are categorized into irrep spaces of 0° in Table 

Harmonic Trap: The largest point symmetry possible in is 0(2), all orthogonal transformations 
of the plane, i.e. reflections through and rotations about the origin. This is the configuration space 
symmetry for the harmonic potential. 1 call this an emergent symmetry because, unlike the cases 


of the asymmetric and symmetric traps, for a harmonic trap the non-interacting configuration space 

symmetry C® cannot be generated by single particle symmetries and particle permutations. 

^ By D4, here I mean the two-dimensional point group, i.e. the dihedral group with four reflection axes that 
is the symmetry group of a square. Coxeter notation for this pure reflection group is BC 2 or [ 4 ]. The same 
symbol D4 is also Schonflies notation for the three-dimensional point group with Coxeter notation [ 4 , 2 ]"*". 
These two groups are isomorphic, but have different geometrical realizations. In the three-dimensional sense, 
the group D4 is an order eight group consisting of only rotations and no reflections. It can be visualized as 
the symmetries of a square parallelepiped with sides two-color checkered by an even number of checks. The 
Schonflies notation for the three-dimensional version of the reflection group D4 is C 4 ,, and it is the symmetry 
of a square parallelepiped with the two square ends painted different colors. 
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Table 1 This table categorizes the energy levels in the two-particle, non-interacting spectrum (J 2 for a general 
symmetric trap into irreps of the non-interacting configuration space symmetry group C® ~ S 2 x 0(1)’^^, into 
equivalence classes of irreps of the kinematic symmetry group K® ~ S 2 x (0(1) x Tt)^^, and into equivalence 
classes of irreps of the interacting symmetry group K 2 ~ S 2 x 0(1) xTt. Non-interacting energy levels are labeled 
by their one-particle compositions. See text for notation for irreps, and see Fig. for the partial ordering of 
these energy levels if the specific values of the one-particle spectrum are unknown. Note that because K® ^ 0° 
and K§ O K 2 , generally the K® irreps are reducible with respect to C% and K 2 . 


Composition 

Degeneracy 

C 2 irreps 

K 2 classes 

K 2 classes 

Lo^j 

1 

Ai 

L+^J 

[2] + 

Loij 

2 

E 

L+-J 

[2]-©[l^l- 

[Pj 

1 

A 2 

L-"J 

[2] + 

L02J 

2 

4li ©Bi 

L+ 1 + 2 J 

[2]+ © [1^^] + 

L12J 

2 

E 

L+-J 

[2]-©[l^l- 

L03J 

2 

E 

L+-J 

[2]-©M- 

L2"J 

1 


L+^J 

[2] + 

L13J 

2 

A 2 0 B 2 

L- 1 - 2 J 

[2]+ © [1^] + 

L04J 

2 


L+ 1 + 2 J 

[2]+ © [1^] + 


The irreducible representations of 0(2) are labeled by m S N and they are one-dimensional for 
m = 0 and two-dimensional for to > 0. They correspond to the polar harmonics exp^imcj)). One can 
think about this as ‘angular momentum’ in configuration space and construct an observable L 12 ~ 
i{dia\ — a 2 a\) out of ladder operators that commutes with Hq. The degeneracy of the total energy total 
energy E = hu!{X + l) is d{E) = X-l-1, so the dimensions of the irreps of 0(2) are insufficient to explain 
the degeneracies of u® foi' fhe harmonic trap. Explaining the total degeneracy requires considering the 
full kinematic symmetry of E[q. 

3.1.2 Two Non-Interacting Particles: Kinematic Symmetries 

The previous section established that the minimal two-particle configuration space symmetry is given 

by 

C° D P 2 X Of 2 (17) 

This expression relates the one-particle configuration space symmetry Ci to the two-particle non¬ 
interacting configuration space symmetry of Hq for both symmetric and asymmetric traps. 
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A similar situation holds for the kinematic symmetry of iJg. For general symmetric and asymmetric 
traps, the minimal kinematic symmetry is 


D P2 >< 


(18) 


The one-particle kinematic symmetry Ki always includes time translation Tt, so now there are two time 
translations, one generated by each particle’s Hamiltonian Uj{t) = exp{—iHjt). Since the particles are 
non-interacting, their clocks are not linked and their time lines are independent. Total time evolution is 
also a symmetry of course, but for non-interacting particles it can be generated by single particle time 
evolutions U{t) = Ui{t)U 2 {t). Note that the exchange operator does not commute with the one-particle 
time-translations; instead one finds U{12)Ui{t) = U2{t)U{12). This means K 2 is not abelian and so its 
faithful irreps are not one-dimensional. 

Asymmetric Traps: For asymmetric traps, Ki is just Tt and the minimal kinematic symmetry group 
is K 0 ^S 2 >< The irreps of K 2 are labeled by the state composition [i/J and the irrep space is 
denoted 1C . The decomposition of the spatial Hilbert space into K 2 irreps is the a direct sum over 
all compositions spaces: 


jC= 0 

[vj Gcri XCTi 


(19) 


Each K, is an energy eigenspace with energy . Unless there are emergent symmetries or accidental 
degeneracies, then each is distinct. 

The irreps of Kg fall into two equivalence classes. Compositions like with shape [v] = 

[2] have one-dimensional irrep spaces spanned by |e@). Compositions \y\ = [a/3J with shape 

[v] = [1^] have two-dimensional representation spaces /CL«/3J spanned by |@2]) and ||]). Note that since 
K§ D CO, irreps of K® may be reducible with respect to C®, for example \_aP\ ~ [2] 0 [12]. 

Symmetric Traps: The inclusion of parity symmetry does not change the irrep structure or change 
the decomposition (19), but now there are five equivalence classes instead of twc0 These equivalence 
classes are: and [—, compositions of two copies of the same state with even parity or odd parity; 

[ 0102 J and [— 1 — 2)5 compositions of two different states both with even parity or odd parity; and 

^ These are not the same five irreps as Cg ~ D4, but the fact that the number of equivalence classes of 
irreps is the same as the number of irreps is valid for any N. 
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[H—J, compositions of an even and odd state. See Table for examples of how compositions spaces 
are sorted into Kq irrep equivalence classes and reduced into Cg irreps for low-energy compositions. 

Harmonic Traps: Two particles in a harmonic trap is the same as a two-dimensional isotropic 
harmonic oscillator, and so K 2 ~ U(2). A representation of u G U(2) is defined by operators U{u) that 
act on the pair of one-particle annihilation operators a = (fii, 02 ) as 

U{u)a.U\u) = ua. (20) 

These transformations leave the non-interacting Hamiltonian Hq = hu!{a}a + 1) invariant. The group 
U(2) is isomorphic to the set of all symplectic, orthogonal transformations of four-dimensional phase 
space. The irreducible representations of U(2) are equivalent to the more familiar SU(2): they are 
finite-dimensional and labeled by an non-negative integer X. This quantum number is the same as the 
total excitation X = ni + n 2 of a pair of oscillator^ The dimension of the irreducible representation 
is d(U(2);X) = A -|- 1. As a consequence, when A > 1 there must be multiple compositions \ v\ with 
the same energy, not just the one- or two-fold degeneracy inherited from the one-particle symmetry 
via the subgroup S 2 x (U(l) x Tt)^^. These degeneracies imply the existence of other operators besides 
those generated by S^r and Ki that commute with Hg. Several inequivalent complete sets of commuting 
operators can be chosen and these correspond to the different coordinate systems in which the two- 
dimensional isotropic harmonic oscillator separates, i.e. cartesian, polar and elliptic m- 

Another kind of emergent two-particle kinematic ‘symmetry’ results from accidental degeneracies. 
The most famous of these are the Pythagorean degeneracies that occur for the infinite square well (see 
for example, |52|1. These are not usually interpreted as symmetries because there is no corresponding 
(linear or non-linear) transformation on configuration space or phase space that induces a unitary 
representation on the whole Hilbert spac^ 

For SU(2) is standard to use j = A/2 as the label. 

® Of course, what one person calls an accidental degeneracy could be an undiscovered symmetry! It seems 
unlikely that after all this time that Pythagorean degeneracies will find a description in terms of configuration 
space or phase space transformations, but there may be other cases of accidental degeneracies waiting to be 
revealed as globally-defined emergent symmetries. 
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Technically one can construct operators which exploit the accidental degeneracy as a symmetry, but 
to do so requires knowledge of the spectrum. For example, the square well state^with (ni, 712 ) = (0, 6), 
(6, 0) and (4,4) span a three-dimensional subspace = /CLoej with energy SOcq. One can define 

a family of operators isomorphic to U(3) that act unitarily on the three-dimensional energy eigenspace 
and act as the identity on the rest of the spatial Hilbert space K. © Those operators would 
realize the accidental degeneracy as a kinematic symmetry group. However, the construction of such 
operators requires prior knowledge of the degeneracy instead of actually explaining how the degeneracy 
arises from the kinematic symmetry of the Hamiltonian and acts trivially on most of /C. It is therefore 
not as useful as a true emergent, global kinematic symmetry. 


3.2 Two Particles: General Two-Body Interactions 


Now add a two-body interaction to the Hamiltonian: 

= Hi + 


( 21 ) 


Only Galilean-invariant two-particle potentials V 12 are considered. The requirement of Galilean invari¬ 
ance can be summarized algebraically in terms of commutation relations: 


[^, V12] = [Qi, V12] = [Q2, V12] = [G(12), 1/12] = 0 
and [A,Vi2] =-[A,I>12]. 


( 22 ) 


The second line of (221 is equivalent to saying that the two-body interaction commutes with the center- 

this implies that the interaction can be written 


of-mass motion. Combined with the first line of 
as V 12 = {'/2\Ri\), where Ri = (Qi — Q 2 )I'J^ is the normalized relative position coordinate. 


3.2.1 Two-Body Matrix Elements 

The condition [[/(12), V 12 ] = 0 also implies the two-particle matrix elements of the interaction {oil3\Vi2\^S) 
have the property 

(a^|Vi2|7<^) = (/3a|Vi2|^7) = VP6\ ■ (23) 

® With the convention that the ground state has n — 0, the energy of infinite square well is = eo(u. -f 1)^- 
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This notation for the matrix elements emphasizes that this amplitude is relevant for the state transitions 
a e-)- 7 and /3 O i5. The one-particle basis can be chosen so that these matrix elements are all real. 
The group P 2 is a symmetry for both and V 12 , so it remains a symmetry of the total interacting 
Hamiltonian . Therefore there are only matrix elements between states carrying the same irreducible 
representation of S 2 , i.e. 

(E33| '1^12 i> = (Effll i> = 0 (24) 


for any states a, /3, 7 and S. 

First order perturbation theory gives the level splitting of the non-interacting states in the limit 


of weak interactions. In terms of the interaction matrix elements for the symmetrized states (111, the 
level splittings are 


(EHI V12 I03) = [/32J -I- 

(II 1^12 II) = V[a^] [/32J - I^[a/3J2, (25) 

where for brevity I denote U[q, 2 j|^^ 2 j = v\^aa\ipi 3 \ and 'C[q,/ 3 j 2 = v\^a 0 \lai 3 \ - This implies the familiar 
result that for two-particle interactions the interference between the direct channel and the exchange 
channel generally shifts the symmetrized state more than than the antisymmetric state. 


3.2.2 Symmetries of Two Interacting Particles 

The minimal non-interacting symmetry P 2 x is partially broken by Vi 2 - The particle exchange 
symmetry P 2 is preserved. The diagonal subgroup of i.e. elements like II 1 II 2 and Uilt)tJi{t), 

still commutes with H^. Therefore the kinematic symmetry of the interacting two particle system K 2 
always contains a subgroup isomorphic to the one-particle kinematic symmetry Ki, so: 

K 2 D P 2 X Ki. (26) 

Asymmetric Trap: The minimal total kinematic symmetry of the interacting system in an asym¬ 
metric trap is K 2 ~ S 2 X T( and the configuration space symmetry is just C 2 ~ S 2 . Both of these groups 
are abelian with only one-dimensional irreps, and so in this minimal case, the spatial Hilbert space 
decomposes into irrep towers K. = 0 Each energy in tT 2 is non-degenerate and associated to 

either the S 2 irrep [ 2 ] or [ 1 ^], unless the interacting system has emergent symmetries. 
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Symmetric Trap: If the trap respects parity then Ki ~ 0(1) x Tj. This group has the same irreps 
as the asymmetric case, but double the number of irrep equivalence classes, so K. is decomposable into 
four towers 

K. = 


The total parity operator 11 commutes with V 12 and parity remains a good quantum number even 
when interactions are turned on. Therefore, in addition to selection rules against transitions between 


states with different exchange symmetries (24), matrix elements of V 12 between two-particle states 
with different parity must be also be zero. This reduces the number of matrix elements required for 
exact diagonalization in a truncated Hilbert. See Table for the reduction of K® irreps into K 2 irreps 
for the case of symmetric traps. 

Note however that the one-particle parities Tli do not commute with V 12 . For symmetric traps 
the configuration space symmetry is reduced from 0 ° ~ D 4 (with order eight and five irreps) to 
only C 2 D 2 (with order four and four one-dimensional irreps). For two particles, the permutation 
operator 17(12) can also be interpreted as relative parity: reflection across the line <71 = 92 reverses 
relative position U[12)Ri = —Ritj[12). The operator 17(12)17 is a reflection across the line qi = —q 2 
that reverses the normalized center-of-mass position R = {Q 1 +Q 2 ) and leaves the relative position 
Ri invariant. 

Harmonic Trap: For the harmonic trap K 2 ~ S 2 x U(l) x T* and the extra U(l) symmetry provides 
an additional good quantum number: the center-of-mass excitation n. The spatial Hilbert space /C is 
decomposable into an infinite number of equivalence classes, one for each value of n, and these each 
further separate into parity towers: 


The total parity of states in is tt = (—1)" and states in 1" have parity tt = —(—1)". The 
selection rules preclude non-zero matrix elements among states in different towers and this makes 

a significant reduction in effort for calculating higher-order terms in a perturbation series or for making 
exact diagonalization in truncated Hilbert spaces. 

In summary, the minimal kinematic symmetry K 2 of the two-particle interacting Hamiltonian 
is P 2 X Ki. This group has only one-dimensional irreps, and so unless there is an accidental or emergent 
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symmetrieQ there are no degeneracies in the interacting spectrum a 2 - The symmetry K 2 is enough 
to completely specify the qualitative features of level splitting for weak interactions. To calculate the 
specific energy shift requires the two-particle interaction matrix elements of the form U|^o, 2 j |^^ 2 j and 
U|^q,^j 2 , but the splitting is universal. 

3.2.3 Two-Body Contact Interactions 


Now, specify V 12 to be the contact interaction, which in the position representation is 


V^ilqi - < 72 !) = gSiqi - 92 ) = -^gS{ri), 


rr 


(27) 


where ri = (gi — q2)l\f2 is the normalized relative position coordinate. This potential satisfies the 


Galilean invariance requirements (22) and therefore the kinematic symmetry group contains at least 


the minimal symmetry P 2 x Ki. The goal of this section is to find results that are trap-independent 
using symmetry methods alone. Note that the case of the contact interaction is analytically solvable 
for two-bodies for any value of g in one-dimensional harmonic trap [5411551156j or for infinite square 
well [HIST]. Finding the energy for general g requires solving a transcendental equation, but the system 
is integrable for both of these traps. 

For the contact interaction, the two-particle matrix elements {aj3\Vi2\')S) = are invariant 

under permutations of the four states a, /3, 7 and S. This is shown by by going to the position 
representation wher^ 

{al3\Vi2\"f5) = J dqldq2^|^UglWfiig2)gSiql - q2)'ip.y{qi)-ipsiq2) 

= g dqipa{q)'4’l3iq)'<pjiq)tpsiq) 


= V 


[alS-yS] • 


(28) 


This ‘state permutation symmetry’ of the contact interaction matrix elements means that in addition 


to zero matrix elements between states in different irreducible representation spaces of S 2 as in (24), 


^ Another example of emergent symmetries is the case of harmonic interactions F^(\/2|Ai|) oc in a 

harmonic trap. This system has interacting kinematic symmetry K 2 ~ (U(l) x because both the center- 

of-mass and relative coordinate act like a one-dimensional harmonic oscillator. 

® Remember that wave functions of the one-particle Hamiltonian can be chosen as real without loss of 


generality. 
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the matrix element between totally antisymmetric states is also necessarily zero 




(29) 


as one shows by inserting n|^Q,2j |^^2j =Vy^p ^-2 =vyop.p 2 ^ in (25). The consequence, as expected, is that the 
fermionic states do not “feel” the contact interaction and remain stationary states of the Hamiltonian 
for all values of the interaction strength g, attractive g < 0 or repulsive 5 > 0. 


3.3 Two Particles: Unitary Limit of Contact Interactions 


In the unitary limit g —> 00, the contact interaction is like a sword through configuration space, severing 
the two halves with a nodal line that no probability current can penetrate. Each sector in configuration 
space acts as a disjoint domain for wave functions. The particles are either in the specific left-to-right 
order qi < q 2 or in the order q 2 < qi- The spectrum in the unitary limit is therefore the same as 
the spectrum of totally-antisymmetric non-interacting states |[|]). There is a two-fold degenerate level 
for every state pair a ^ j3 (see the right side of Fig. [^. As in the non-interacting case, only a 
partial ordering of the spectrum can be determined without specific knowledge of cti. 

Define the ‘snippet’ basis wave functions [inillll for each pair a yf /3 and each order qi < qj by 


(q| |;{i 2 }) 

(q| |;{2i}) 


V^(q||) = </'a(9l)</'/3(<72) - (t>a{q 2 )<l) 0 {qi) for qi < q 2 

< 

0 for q 2 < Qi 

\ 

/ 

-V^(q||) = 4'a{q2)4'l3{qi) - (I>a{qi)(t>l3{q2) for q 2 < qi 

< 

0 for qi < q 2 


(30) 


The two vectors ||];{12}) and ||]; {21}) form a basis for the energy eigenspaces of the unitary-limit 


Hamiltonian with energy Cq, + C/?. From (30) they transform under particle exchange like 


U(12)||;{12}> = ||;{21}). 


(31) 


The states that are symmetric and antisymmetric under particle exchange are 

li^Ea) = (li; {12}) -b ||; {21})) 

li;S) ^ ^ (Is {^^{) “ Is {^^{)) • 


( 32 ) 
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The state |[|];[|j) is in fact just the original non-interacting, antisymmetric basis vector |||j^ and the 
state ||]; is its symmetrized version. Although these two states have the same position probability 
density | (q| |]:^) p = | (q| [|]; [|j) P, they will have different momentum distributions because of the 
cusp in (q| gEs). 

When the trap is parity symmetric and the one-particle states |a) and |/3) have parities and Tip, 


then one infers from (301 that 


fl ||;[I11]) = -TTaTT/J ||;E0) 


(33) 


In other words, at unitarity the symmetric state always has opposite parity to the antisymmetric state 
from which it is constructed. 


3.3.1 Two Particles: Near Unitary Limit 

What about the not-quite-unitary limit? Consider a weak perturbation T of IT^ that mimics the 
effect of not quite having an infinite barrier. Such operator would allow a little tunneling between 
the two sectors of configuration space, and should decrease the energy of states like |[|];[is) due to 
the less dramatic cusp at the nodal line qi = ( 72 - It also must have zero matrix elements between 
antisymmetric states because those states do not feel the contact interaction. An operator T that 
satisfies those requirements has uniform matrix elements in the snippet basis: 

(i { 12 }| f II; { 21 }> = (|; { 21 }| f ||; { 12 }> = (|; { 12 }| f ||; { 12 }) = (|; { 21 }| f ||; { 21 }) = -t. 

(34) 

For a specific trap with known energy eigenstates, the small positive constant t can be calculated [29] 
from the wave function (q| 1 ) as 

2 

(35) 

91=92 

As expected, the eigenstates of +T are also the symmetrized |[|];[I@]) and antisymmetrized |[|];|ij) 
states, now with eigenvalues ea + ep — 2t and £« -I- e/ 3 , respectively. 

As another application of this section, by combining the results for weak splitting from cr^ and 
for not-quite-unitary splitting from cr^, a one-to-one adiabatic mapping from non-interacting states 


r*+oo 


t = 


9 J-o 


dq2 


d (q| I) 


dqi 
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to unitary states can be determined, the simplest case of the famous Fermi-Bose mapping |53) . The 
non-interacting symmetric state | m ) is mapped to |[|];on2]^ where fi = a + 1 while the non-interacting 
symmetric state | 0 a) with a < /3 is mapped to |[|;C[I1]) where 7 = /3 -|- 1. This amounts to adding one 
nodal line to each of the symmetric states at the location of the contact interaction. Of course the 
antisymmetric states |||]) = |[|];[ij) are unchanged under the adiabatic mapping because they already 
align with the nodal line at qi = 92 - 


3.4 Spin and Symmetrization for Two Particles 

Before moving on from two particles, let us finally consider the incorporation of identical particle 
symmetrization and spin degrees of freedom and state some well-known results. When there are J > 1 
spin components accessible, then the total total Hilbert space % = JC® S \s the direct product of the 
spatial Hilbert space K ~ and the spin Hilbert space S ~ . The P 2 symmetry implies that 

the total Hilbert space can be decomposed into symmetrized subspaces 

If there are no spin degrees of freedom, then S is one-dimensional and the symmetric subspace of the 
spatial Hilbert space ^ is available for population by identical bosons and the antisymmetric 
subspace ~ by fermions. If there are spin degrees of freedom, then S can also be decomposed 
into symmetric and antisymmetric subspaces and 1. For example, for two particles with spin 
1/2, the triplet states are in and the singlet state is in The spin Hilbert spaces and spatial 
Hilbert spaces are then combined as 

= (x:[ 21©5P1)© and 

•^[P] ^ ^ ^[P] j ^ ^ ^[2] j ^ (-gg) 

Note that the total spin operator S = (Si 4 - 82 )^ and total spin component operator = Siz + S2Z 
commute with the permutation operator 1 /( 12 ), as well as with all one-particle and two-particle spatial 
observables. Therefore, total spin and spin component are good quantum numbers for the symmetrized 
states for any interaction as long as the trap is spin-independent. Symmetrization induces correlations 
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between spin states and energy, for example making the lowest energy state two spin- 1/2 fermions only 
accessible to the singlet combination. 


4 Three Particles 


The kinematic symmetry of the Hamiltonian (whether interacting or non-interacting) includes 
particle permutation symmetry P 3 ~ S 3 . Unlike S 2 , the group S 3 is not abelian, and so now the 
irreducible representations of particle exchange symmetry are more complicated. One implication is 
that the spatial Hilbert space JC, the spin Hilbert space S and the total Hilbert space "H can each be 
broken into subspaces with the three types of symmetry that three particle states can have, e.g. for 
the spatial Hilbert space 

/C = (37) 


Each of these subspaces in (37) is a tower of irrep spaces of S 3 . If S 3 is the only symmetry of 77^, 
then each copy of each S 3 irrep would have a distinct energy. 

One complication of S 3 compared to S 2 is that the elements of S 3 are not all their own inverses, 
and so now we have to be a little more careful about representations. I choose the convention that a 
particle permutation p € P 3 acts on the coordinates like 


Qip){qi,q2,q3) = {qpi, qp2, Qpi), 


(38) 


where p is expressed in permutation notation {piP 2 P 3 }- Using ^^(q) = (qltf'), the induced representation 
on wave functions is 

U{p)^{q) = 1'{0{p~^)q). 


This implies that particle permutations are represented on the particle basis |nin 2 n 3 ) similarly to 


eq. (38): 


U{p)\nin2n3) = \np^np^np^). 


(39) 


This convention is opposite to the convention used in |41j . which otherwise (especially chapters 3 and 
4) provides an excellent reference for the methods used in this section and in the sequel. 

Similar to the previous section, the following subsections treat the cases of non-interacting, inter¬ 
acting, and contact interactions in the unitary limit. The usefulness of state permutation symmetry 
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Fig. 4 Each sequence of three numbers determines a composition [v\. The superscript denotes the number 
of sequences with that composition, or equivalently the dimension dlu\ of Unless the specific values of 

ti € (71 are known, only the partial ordering of composition energies given by the arrows is dehned. The 
boxed sequences have compositions with the shape [1^]. Their composition subspaces carry the regular 

representation of S 3 and have state permutation symmetry ~ S 3 . There is a one-component bosonic 

state in every composition space K. , but there are only one-component fermionic states in spaces with the 
boxed compositions. Note that if the sequence (0,1, 2) is added element-wise to each of the original sequences 
(e.g. (0, 0, 2) -I- (0,1, 2) = (0,1,4)) the chart will have the same form, giving a one-to-one mapping from bosonic 
states to fermionic states, as in Fig. 


becomes more evident as the limits of particle permutation symmetry become more acute, and in the 
unitary limit of the contact interaction a new kinematic symmetry emerges called ordering permutation 
symmetry. 


4.1 Three Particles: Non-Interacting 

The spectrum of three non-interacting particles dg is constructed from all possible compositions of the 
single particle energies. As with two particles, the spectrum dg can be only partially ordered without 
specific knowledge of di (see Fig. |^. If there are no emergent or accidental symmetries, then there 
are three kinds of energy levels: singly-degenerate levels derived from compositions of identical states 
like [q:^J , three-fold degenerate levels from compositions of two different states like [a^/3J, and six-fold 
degenerate levels from compositions of three different states like [a/JyJ. Since S3 has only one- and 
two-dimensional irreps, it is clear that S3 symmetry alone cannot explain the degeneracies of dg. 
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The configuration space and kinematic symmetries of three non-interacting particles inherit the 
following subgroups by their tensor product construction: 

D P 3 >< Cf 3 

D P 3 X (40) 

Indeed the group P 3 x has irreps that are one, three and six dimensional, as shown below. 
These irreps are labeled by the the three energies in the composition, i.e. the three characters of the 
time translation subgroup T^^^. Therefore, the minimal kinematic symmetry is sufficient to explain 
the degeneracy of the non-interacting energy levels in Ug unless there are emergent symmetries or 
accidental degeneracies. 

4 .. 1.1 Three Non-Interaeting Particles: Configuration Space Symmetries 

Asymmetric Trap: The minimal configuration space symmetry occurs when Ci ~ Zi and C 3 ~ S 3 . 
The equivalent point group in three dimensions has Schdnflies notation Cg^,. This is the symmetry of a 
triangular prism with distinguishable ends. The permutations p G Pg of three particles are realized in 
by an orthogonal matrix 0{p) S 0(3). Each of the three two-cycles 0{ij) is a reflection across the 
plane defined by Qi = qj and the two three-cycles 0(123) and 0(132) are rotations by ±27r/3 about 
the line 91 = <72 = 93 - See Fig. for some examples of equipotential surfaces for three non-interacting 
particles. 

Symmetric Trap: When the trap is parity symmetric, each particles’ parity operator remains a 
symmetry of the system. Then Ci ~ 0(1) and the configuration space symmetry is 

C°~S3><0(1)^3_0^, (41) 

where Oh is the full cubic symmetry in three-dimensions with order 48 and ten irreducible repre¬ 
sentations [SS], five with even parity and five with odd. See Table for a categorization of low-level 
non-interacting three particle states using the standard notation [50] for the Oh irreps Aig, Eg, 
etc. 

Harmonic Trap: For a harmonic trap, the configuration space symmetry is maximal: Cg ~ 0(3). 
As a result, another good basis of energy eigenstates is provided by the quantum numbers {X, A, m}. 
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Fig. 5 These figures depict a single equipotential for three particles in a one-dimensional trap with (a) V(q) = 
(b) V(q) = kl; (c) V(g) = q^ (harmonic); (d) V(q) = q^°; (e) V(q) = kl for 5 < 0 and V(g) = q^° 
for q > 0. The three transparent planes are the surfaces in configuration space where two particles coincide, 
i.e. qi = q 2 , ?2 = 93 and q\ = gs. All subfigures have particle exchange symmetry P 3 : two-cycles like (12) 
are reflections across the corresponding coincidence plane and three-cycles like (123) are rotations by 27r/3 
about the axis where all three planes intersect. Potentials (a)-(d) are also parity symmetric and have the full 
octahedral symmetry Oh- Note that only (c) the harmonic potential (which has the maximal point symmetry 
0(3)) is symmetric under relative parity inversion, corresponding to rotation of tt about the bold axis where 
the three coincidence planes intersect. 

where X is total excitation, A € N labels the 0(3) irrep and is like orbital angular momentum, and 
m S {—A,..., A} labels an orthogonal basis within the irrep. Alternatively, the subgroup 0(2) can be 
used to decompose the spatial Hilbert space into cylindrical harmonics and a center-of-mass quantum 
number. The cylindrical basis has proven particularly useful for exact diagonalization when interactions 
are included [53J [5SJ |35] . 

4 .. 1.2 Three Non-Interaeting Particles: Kinematic Symmetries 

As for the two-particle case, the irreps of the minimal kinematic symmetry Kg ^ P3 >< are labeled 
by compositions \ v\ . In the simplest case of an asymmetric trap, these irreps fall into three equivalence 
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Table 2 This table categorizes the energy levels in the three-particle, non-interacting spectrum cr® for a general 
symmetric trap into irreps of the non-interacting configuration space symmetry group C3 ~ S3 ix 0(1)^^ ~ Oh, 
the kinematic symmetry group K3 ~ S3 ix (0(1) x Tt)^®, and the interacting symmetry group K3 ~ S3 x 0(1) x 
Tt. Non-interacting energy levels are labeled by their one-particle compositions. The notation for Oh irreps 
is standard, c.f. [58], and the notation for the irreps of K 3 is the composition with subscripts denoting each 
state’s parity. See Fig. for the partial ordering of these energy levels if the specific values of the one-particle 
spectrum are unknown. One C 3 irrep does not appear in this table; the lowest energy composition that carries 
a copy of the irrep Eu is [1^3J —>■ A 2 U © Eu- Two K 3 irreps do not appear in this table; the lowest energy 
compositions that carry those irreps are [024J —>• [+1 +2 + 3 J and [135J —>■ [—1 —2 — 3 J. 


Composition 

Degeneracy 

Co irreps 

Kq irreps 

irreps 

Lo^j 

1 

Alg 

L+^J 

[3]+ 

Lo"iJ 

3 

Tiu 

L+^-J 

[3]-©[21]- 

Loi^j 

3 

T2g 

L+-^J 

[3]+ © [21]+ 

L0^2J 

3 

Alg 0 Eg 

L+?+ 2 j 

[3]+ © [21]+ 

Li^j 

1 

A2u 

L-^J 

[l^]- 

L012J 

6 

Tlu © T2u 

L+i +2 -J 

[3]- ©2[21]- ©[l^]- 

[0^3] 

3 

Tiu 

L+"-J 

[3]-©[21]- 

L1"2J 

3 

T2g 

L+-^J 

]3]+ © [21] + 

[02^) 

3 

Alg 0 Eg 

L+?+ 2 j 

[3]+ © [21] + 

L013J 

6 

Tig © T2g 

L+ -1-2J 

[3]+ © 2[21]+ © [1^] + 

[0^4] 

3 

Alg © Eg 

L+H2J 

[3]+ © [21] + 


classes, depending on the shape of the composition. One basis for irreps spaces AiiLO is the particle 
basis of all sequences in the composition, i.e. 

J = span{|aaa)} 

= span{|aa/?), |a/3a), |/3aa)} 

^La/ 37 J = span{|a/ 37 ), |a 7 / 3 ), \Paj), l/S-fa),\jal3),\jl3a)} . (42) 


In the basis (42), the subgroup has been diagonalized. A complete set of commuting observables for 
this basis are the generators While this basis is natural, it is not optimized for the tasks 

of perturbation theory and exact diagonalization when interactions are incorporated, or for analyzing 
the unitary limit of contact interactions, or for symmetrization of identical fermions or bosons with 
or without spin. For all these tasks, a basis optimized for the S 3 subgroup works better. The rest of 
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this subsubsection describes a particular choice for that basis and the corresponding complete set of 
commuting operators. 

Asymmetric Trap: Continuing with the asymmetric trap, the first step is to reduce the Kg irreps 
into S 3 irreps labeled by [/x]. For each of the three classes of Kg irrep spaces, the reduction looks 

like 




^La"/3J ^ 
^La/37j = jC 





(43) 


In each reduction, the irrep space is reduced into subspace labeled by a semi-standard 

Weyl tableaux W. To make a Weyl tableaux, a Ferrers diagrams is filled with the state labels in the 
composition. These labels must stay the same or increase to the right and must increase to the bottom 
(assume a < /3 < 7 ). There is only one way to fill a Weyl tableau out of the composition [a^J, two 
ways for the composition [a^j3\, and four ways for \ aj3^\. Each subspace C /CLC is isomorphic to 
the S 3 irrep space with the shape \W] of the tableau W. Note that the dimensions of the irreps 

tally as they should in ( [4^ , for example in the last line ^[a/SyJ = d[3] -I- 2d[21] -|- d[l^] = 6. The proof 
of this reduction and the extension for N particles is found in the sequel article. The key observation 
is that K^ irrep spaces are permutation modules [33] of Sat characterized by the shape of the 
composition [v\. I now consider each equivalence class in turn. 

First, irreps spaces like carry the trivial, totally symmetric representation and have 

the single basis vector |iEE]) = \aaa). 

The second equivalence class of spaces K are those like /C . This representation (also called 
the defining representation) is reducible into a totally symmetric sector ~ and a sector 

with mixed symmetry ^ Following [IT], The three basis vectors of can be chosen 

as 

| f^) = (|aa/3) + \al3a) + \(iaa)) 

v3 

iPi®) = ^ (2|aa/3) - |a/3a) - |/3aa)) 
liPiP) = {WPa) — \l3aa)). 


(44) 
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The second basis label for the [21] irrep are the only two Young tableaux with the same shape as the 
Weyl tableau in the first basis label. When filling the Ferrers diagrams for a standard Young tableau, 
the particle numbers must increase to the right and to the bottom. The state | m^ ) does not need a 
Young tableau label because there is only one Young tableau [iMai with shape [3], so it is understood. 


These states (44) are the simultaneous eigenvectors of and two other operators: 


^ ^^^ 2 ) + [7(23) + 17(12) and = 17(12). 


(45) 


The notation here seems a little chunky for three particles, but it can be generalized to N particles 
and to other kinds of permutations, so I have adopted it. The operators are class operators over 
the all c-cycles of permutations of the composition \ v\. Specifically, is the sum of all two-cycles 


' I 12 I 

in P 3 = Pli 23 J and C 2 is the only two-cycle in P 2 = P[i 2 j • The basis vectors (44) are eigenvectors 


of the operators in (45): 


|k^) = 3 |msi), 

1 

f™) = , 

IIPiF) = 0 , 


\\mm/ ~ IterteF/ ’ 

IPP) = 0 , 

ci'“ 

— _ |mp^\ 


The two states |||P|P) and |||P|i|®) are therefore distinguished by whether particles 1 and 2 are 
symmetric or antisymmetric under exchang^ Therefore the set are commuting 

operators that diagonalize composition spaces like /CL“ with respect to the S 3 subgroup of Kg. 

The third type of composition space K is six-dimensional and carries the regular representation 
of S 3 , i.e. each representation appears as many times as its dimension. That means that the space 
appears twice and the set of commuting operators {H^, } does not completely diagonalize 

K. . 

One way to solve this problem is to introduce state permutation group ^[jyj. The state permutation 
group ^[q,/ 37 J is isomorphic to S 3 contains the six elements e, (a/3), (Pj), (ay), (a/ 37 ), {ajP). 
State permutations exchange state labels, not particle labels. Note the following comparisons that show 
® This is equivalent to diagonalizing the S 3 by the canonical subgroup chain S 3 D S 2 and up to a phase the 


basis (441 is equivalent to the Yamanouchi basis |41l I50| . 
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they are distinct operations on K, : 

U{12)\aj3^) = I/Say) = J 7 (a/ 3 )|a/ 37 ) = I/Say), and 
U{V2Z)\aj3^) = l/Sya) = J 7 (a/ 37 )|Q!/ 37 ) = \l5^a), but 
^ 7 ( 12 ) 1 / 370 ?) = |7/3a) ^ IJ{al3)\f5^a) = |a7/3), and 

f/(123)|/3a7) = |a 7 ^) ^ U{af3^)\l3a'i) = |7/3a). (46) 

The elements p G P 3 and p € 'P \_a0-y\ commute with each other on K, , and there is a one-to-one 
map between the groups established by their action on the intrinsic state jo/Sy) |41j . 

State permutations are not a symmetries of Hg. Exchanging two states changes in a particle basis 
vector usually changes the energy of the state. However, the state permutations '^iv\ of a composition 
\v\ are defined so that they leave the spaces invariant. For compositions like and [a^/3J, 
no state in those composition appears the same number of times and so state permutations are trivial 
~ ~ Zi. State permutations fP[c(/ 37 j distinguish the degenerate [21] irreps of S 3 in /CL“/ 37 j 

by introducing the operator 

= U{aP). (47) 

Then, choosing the phase convention of [41], the following basis vectors are simultaneous eigenvectors 
of the set operators which do not commute on all of /C, but do commute 

on /CL“/ 57 J. 

jEES) = ^ (|a/?7) + I/Jay) + l7/?a) + |a7/3) + [ 70 ^) + l^7a)) 

(2|q;/37) + 2|/3q! 7) - |7/3a) - \a'yj3) - jya/S) - j/Sya)) 

\ (-|7/3a) + lay/?) - lya/?) -f \ha)) 
liPiF) = \ (“lT/?a) + lay/?) + |ya/3) - |/3ya)) 

|iPi®> = (2|a/3y) - 2|^ay) -P |y/3a) -7 |ay/3) - |ya/3) - |/3ya)) 

^ (|a/?y) - |/?ay) - |y/3a) - |ay/3) -7 |ya/3) -7 \l3-fa)). (48) 

As before, the the shape of the Weyl tableaux W denote an irrep of S 3 and the Young tableaux Y G [W] 
identify a basis for the irrep spaces. Particle permutations p G P 3 mix states with the same W and 
different Y’s, and state permutations p G ^iap-y\ mix states with different W’s and the same Y. In 
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Chapter 3 of |41] . an explicit method for constructing a complete set of commuting operators that 
correspond to the W and Y is described for the general case of N particles. 

Symmetric Traps: For symmetric traps, irreps of Kg are still labeled by compositions, but including 
parity information there are now ten equivalence classes. There are two equivalence classes [+^J and 
[—with the shape [3]; four classes [+ 1 + 2 J, L+^~Ji L~i~ 2 j with the shape [ 21 ]; and four 

classes [+1 +2 + 3 J, [+1 +2 —J, [H —1 — 2 J and [—1 —2 — 3 J with the shape [l^j. These reduce into S 3 
irreps the same way ( |43[ ) as the asymmetric trap. See Tab. [^for examples of low-energy compositions. 

Harmonic Traps: For the case of the harmonic trap, a state with total energy E = huj(X + 3/2) 
has a degeneracy d(U(3);X) = [X -\- 1){X + 2)/2. These additional coincidences in cTg are explained 
by the emergent U(3) kinematic symmetry |601161j and there are multiple inequivalent complete sets 
of commuting observables that separate the spatial Hilbert space. 

4.2 Incorporating Spin and Symmetrization 

For non-interacting spinless particles, the spectrum and degeneracy is now effectively solved. Since 
~ spinless distinguishable particles can populate every energy level in Ug. Identical spinless 
bosons are restricted to the sector ^ where is the tower composed of the single, spatially- 
symmetric state that exists in every composition space JC . Spinless fermions are restricted to the 
sector ^ spatially antisymmetric states exist only in composition spaces like . There 

is a one-to-one mapping between the set of all composition spaces and the set of composition spaces 
with shape [1^]; see Fig. 

If there are internal components, there are two methods for symmetrization often used. One method 
fixes the spin components of specific particles, e.g. “particle I and particle 2 are spin up and particle 
3 is spin down.” Then the spatial wave functions are symmetrized within particles with the same spin 
components. Examples of this approach include [nilllllilllSlIlSlIMlIMlEZlIMlISn]- The alternate 
method pursued here is the combined, simultaneous symmetrization of spin and spatial states, cf. 
[5511^. Following that approach, the spin Hilbert space S ~ can be decomposed into subspaces 
with definite symmetry 


5 = 
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Then the bosonic sector of % is 


■^[3] = (/c[3l (g) 

[3] 


(49) 


where jp] means the one-dimensional, symmetric subspace of ©/Cl^^l. Bases for this subspace can 
be explicitly constructed using Clebsch-Gordan coefficients for S 3 , c.f. [HIISO]- A similar expression for 
the fermionic sector is 


^[1^1 = © 5[31) © (/CI211 © 51211) © (/Cl^l © 5[i'l). 


(50) 


Note that 5[i^l is non-empty only if J > 3; there must be at least three spin components in order for 
the spin state to carry the required antisymmetry to balance a totally symmetric spatial state. 

Explicit state construction requires additional algebra, but counting degeneracies does not. As an 
example, consider three fermionic spin-1/2 particles with state labels f and ], corresponding to the 
eigenvectors of the z-component of each particle’s spin Si^^- The four states with total spin s = 3/2 
correspond to the totally symmetric spin vectors that span 5131 that are labeled by the Weyl tableaux 
with shape [3]: 

[ pTFiTn ) , Iefeted) , Ieefeb) , Iqzeced) . 

In other words, for three spin-1/2 particles the space 5131 carries one copy of the four dimensional SU(2) 
representation The four states with s = 1/2 that span 512^1 can be chosen as simultaneous 

eigenvectors of 52 , 5^, and 17(12): 

iPFiBpL^x |pnti[i]i]\ i FFnini^ x iMpcsiX 

The first two have 5z-eigenvalue +1/2, the second two —1/2. The space 51231 carries two copies of 
the SU(2) irrep These two copies are distinguished by how they transform under 17(12), as 

indicated by the Young tableau EP or P. 

Combining these results with the spatial symmetries, for every composition space in the same class 
as ^ there are four s = 3/2 fermionic states from the product /C™©5l3l and four s = 1/2 fermionic 

states from the reduction of the product (/C^ © /C^) © 512^1. For every subspace like JC there are 
two s = 1/2 fermionic states in the reduction of /C®^© 512^1. Three spin-1/2 fermions cannot populate 
energy levels like /Cl“ 1 because the spin Hilbert space cannot ‘carry’ enough asymmetry to balance 
the symmetric state. 
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4.3 Three Particles: General Interactions 

Now we add the pairwise interactions 

= Hi + 1/3^ (51) 

where 

1/3 = l/l2 + 1423 + 1^31 • (52) 

In principle, there could also be an intrinsic three-body interaction that satisfies Galilean invariance 
and cluster decomposability, but here I only consider three-body interactions that result from pairwise 
interactions. 

By construction, the operator has P 3 symmetry. It also inherits all the consequences of Galilean 
invariance from the two-particle interaction. Specifically, it commutes with the (normalized) center- 
of-mass position operator and total momentum operator, as well as the total parity operator H = 
/ 71 / 72 /I 3 . It also commutes with the relative parity operator iT^, which is an operator that cannot be 
generated from one-particle symmetries. For two particles, the relative parity operator acts identically 
to the permutation operator {7(12), but for three particles relative parity is not an element of the 
particle permutation group P3 either. To see this, choose the particular set of normalized Jacobi 
coordinates 

R = +Q 2 + Qs) 

Ri = ~ ^ 2 ) 

R2 = ^iQi+Q2-2Q3). (53) 

V6 

Total parity inverts all three coordinates 

HR = — RH^ HR\ = — RiH^ HR 2 — — R 2 H 
whereas relative parity commutes with R and inverts only the relative positions 

tlrpR = RHrp^ Hj-Ri = — R\Hj~^ Hj-R2 = — R2Hr- 

In three-particle configuration space, H^ is realized as a rotation by tt around the line qi = q 2 = Qs- 



41 


4-.3.1 Three Interacting Particles: Configuration Space Symmetry 

Putting this together, the operator has the configuration space symmetry group isomorphic to 
S 3 X 0(1) X (0(1) X T^j). The first factor is due to particle permutation symmetry, the second is 
total parity, and the third is translations and reflections along the center-of-mass axis. The trap will 
certainly break the translational symmetry of V^, and may also break other symmetries. 

Asymmetric Trap: For an general asymmetric trap, the configuration space symmetry remaining 
after two-body interactions are included is only C 3 ^ S 3 ; all parity symmetries are lost. This point 
group is the same as C 3 for an asymmetric trap. 

Symmetric Trap: For a general symmetric trap, the total parity remains a good symmetry so 
C 3 ^ S 3 X 0(1). The Schonflies notation for this three-dimensional point group with order twelve is 
Dsd, Coxeter notation is [[3]], and this is the symmetry of regular hexagonal prism with even-checkered 
sides. This is a subgroup of Cg ^ Oh- 

Harmonic Trap: For a harmonic trap, relative parity provides another independent quantum num- 
beip^ The configuration space symmetry C 3 ~ S 3 x 0(1) x 0(1) is isomorphic to the three-dimensional 
point group and Coxeter notation is [[3], 2]. This is the symmetries of a regular hexagonal prism 
and it is not a subgroup of Oh- 


4-3.2 Three Interacting Particles: Kinematic Symmetry 

Unless there are emergent symmetries, the kinematic symmetry of is Kg ^ Pg x Kg. The irreducible 
representations of this symmetry groups have the same dimensions as the irreducible representations of 
S 3 , independent of the particular one-particle symmetries Kg. Every energy level if G ug is associated 
to an S 3 irrep and the spatial Hilbert space is decomposable into singly-degenerate levels for the totally 
symmetric irrep tower and totally antisymmetric irrep tower and doubly-degenerate energy 
levels for the irrep tower with mixed symmetry Any other degeneracy pattern in erg signals an 

emergent symmetry or accidental degeneracy. 

Relative parity is also a good quantum number for uniform and linear traps because for any quadratic trap 
the center-of-mass and relative coordinates separate. 
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One application of this decomposition is facilitating exact diagonalization in the non-interacting 
basis. Only basis vectors from the same S 3 irreps with the same Young tableau will have non-zero 
interaction matrix elements: 

{WY\V^\W'Y') = {W\\V^\\W')5yy'■ (54) 

This allows the number of basis vectors needed to achieve a certain accuracy to be reduced. For 
example, consider the 56 states that are depicted in Fig. Of those, sixteen states are in and 
only four states are in The remaining 36 states are in but since the interaction operator 

acts like the identity within the irrep, only eighteen states are necessary for exact diagonalization. 
If Ki includes parity symmetry, this provides an additional quantum number tt. Although parity 
does not change the degeneracy of the spectrum 0 - 3 , it can be used to further decompose the spatial 
Hilbert space into two independent irrep towers for each S3 irrep, one for each parity: 

K. = 0 (55) 

Continuing the same example based on the states depicted in Fig. the number of states needed to 
do exact diagonalization in each of these sectors is further reduced to 7, 9, 7, 11, 1, and 3, respectively. 

For harmonic traps, a further reduction is possible. The factor of U(l) in Ki provides an additional 
quantum number conserved by the interactions: the center-of-mass excitation n. The irreps can be 
labeled by To calculate the spectrum 0-3 requires even fewer states because only the case n = 0 
needs to be calculated, e.g. there are only three basis states in Fig. that contribute to the ground 
state in irrep [3 ]q . However, the double tableaux basis \ WY) is no longer the best basis for calculation. 
Instead, observables that exploit separability in cylindrical coordinates on oriented along the center- 
of-mass axis provide a more useful basis 

4-3.3 Three Particles: Weak Interactions 

For weak interactions, the energy of the single state in each space like shifts and the energy levels 

in the spaces like and split and shift. Unless the trap or interaction have additional 

symmetries (emergent or accidental), the degeneracy of the splitting is determined by the symmetry. 
Since Kg D K3, the irreps of Kg are generally reducible with respect to K3. In Table the reduction 
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of composition spaces into K 3 is given. This reduction has a similar form as the reduction of Kg irreps 


by S 3 described by (43). For '' and 0\ the reductions and splittings are identical, and no 
additional information about the nature of the two-body interaction or trap is required. However, for 
compositions like [a/SyJ, the manner in which two copies of split requires specific knowledge of 

the two-body matrix elements. 

Let us make this explicit for each of the three types of composition spaces. In the particle basis, 
the matrix elements of can be expressed in terms of the two-particle matrix elements: 


(a/37|T43|Cr70)=UL„cj [/3?7j ^70 4” ^[/377j [ 70 J 4“ ^[aCJ l7^J 


(56) 


Applying this, one-dimensional spaces like K\'^ ^ experience an energy shift 


(l ° l ° l ° l| |l ° l ° l ° l) = 5viaa][aa] = 


(57) 


The factor of three represents the fact that for this totally symmetric state the pairwise interactions 
of the three particles interfere constructively. 


In three-dimensional spaces like A , the only non-zero matrix elements in the basis (44) are 




(l ° l ° | g || r |l °7 | g |) = W[q,2j2 -I- 2U[q,2J yp2^ + 

(ppi ipr> = (i^ii 


= +2 ULq 2J|^^2J -U|^„^J2. 


(58) 


The non-degenerate symmetric level always experiences a greater shift than the doubly degenerate 
mixed-symmetry level. Additionally, for contact interactions there is state permutation symmetry of 
the two-body matrix elements. Then we have yp 2 y = = U|^c, 2 ^ 2 j and these relations simplify 

further. 

Can anything be considered ‘universal’ for weak perturbations of composition spaces like IC^ 

Yes: the relation between the particle basis and the symmetrized perturbation eigenbasis, the larger 
energy shift of symmetric state compared to the partially symmetric state, and the algebraic expression 
for the energy shift in terms of two-particle interaction matrix elements are all universal features of 
/Cl“ P\ spaces. Although specific numerical values for the weak interaction energy shift depend on the 
specific two-body interactions, those properties do not. 
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For the composition spaces, one of these universal features is lost because such subspaces 

are not simply reducible by S 3 . The reduction of Kg irreps into Kg irreps is still sufficient to determine 
the level shifts in the totally symmetric and totally antisymmetric spaces. The matrix elements of 
in terms of the two-body matrix elements are: 


( 02 q| | 02 h) 





= L/3"J + '^La/3J^ + + ^L/37j^ + h^l + 

= 'f^KJL/3"J -'^La/3J" +^L/3"JL7^J “^L/37J" +^L“"JL7"J “^L«7J"- 


(59) 


These are the first-order energy shifts, and the first-order energy eigenstates remain |05S) and 
as before, in the case of contact interaction the fermionic state will feel no energy shift at this (or any 
order) of perturbation theory. One way to think about this is destructive interference between the direct 
and exchange channels. Symmetry analysis gives another interpretation: because the interaction is 
symmetric under state permutation and the fermionic state is antisymmetric under state permutation, 
matrix elements between fermionic states are all identically zero. 

However, for the two copies of the Sg irrep space in the expansion of , there are are 

matrix elements of between states in and with the same Young tableau. The reduced 

matrix elements are 


And 


3||f3||03\ = 


\14J ll 1114 


all 


1 1 

+^L/3TL7=J - 2^1/371" + [7^1 “ 2^L“7j^ 

1 1 

L/32J - + ^L/3"J 17=^1 + 2'^L/37j" + 17^1 + 2^La7j' 


V3 

= -^L/37jO- 


(60) 


Note that the state permutation symmetry of is broken by the interaction. The state permu¬ 
tation class operator = U{aP) no longer provides a good quantum number and Weyl tableaux 

are not good basis labels for interacting state^^ To first order, the energy shifts for these levels are 
found by diagonalizing in the [ 21 ] sector to find 


'f± = [/ 32 J -I- Vyp2y y^2y + U|^q,2j y^2y 

^\J “ ^La/3J = 3'L/37j" + (^L/37j")^ “ 2 V 2 . (61) 


Also note that the first two equalities of (601 are not equivalent under exchange of /3 and 7 because in the 


choice of basis (48 1 , the state permutation subgroup generated by U{aj5) was diagonalized. 
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Fig. 6 This figure depicts the level splitting under weak contact interactions for a harmonic trap (left) and 
a hard wall trap (right). The energy scale for the two traps has been chosen so they have the same difference 
AE = F|^oi 2 J ~ ^^[ 03 ] between the ground state with composition [O^J and the energy level with composition 
[012J. The strength of the contact interaction is g — AE/30. The thickness of lines depicts the degeneracy of 
the energy level, or equivalently the dimension of the corresponding S 3 irrep. For the highest non-interacting 
energy level, only the level corresponding to the unperturbed totally antisymmetric spatial state is depicted 
and the rest of the levels are cut off. 


The magnitudes of these shifts are intermediate between the shifts in the totally symmetric and totally 
antisymmetric sectors. The corresponding eigenvectors also depend on the two-body matrix elements, 
but the algebraic form of the eigenvalues and eigenvectors in terms of two-body matrix elements do 
not. 


To summarize, for weak interactions, a six-fold degenerate composition subspace like gen¬ 

erally breaks into four levels. The biggest energy shift takes place for the totally symmetric state 
| ragra ) and the smallest energy shift for the totally antisymmetric state See Fig. for an example 
comparing the level splitting of the lowest few energy levels of a harmonic well and a hard wall well 
under weak contact interactions. In contrast, the states |[|]®F) and |||Pl^) mix under the interaction. 
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Unlike the other composition subspaces K, and K, or the totally symmetric and antisymmetric 
sectors of this composition subspace K , specific knowledge of the matrix elements is required 


to determine how the states with mixed symmetry split. The algebraic expression (611 for the level 
splitting of the two partially symmetric subspaces (which are relevant for bosons and fermions with 
J > 2) is universal, but only in the weakest possible sense. In the sequel, it is hypothesized that for five 
particles and more, even this weakest kind of algebraic universality is broken because the diagonalizing 
the perturbation requires solving a quintic equation. 


4.4 Three Particles: Two-Body Contact Interaction at Unitarity 

To visualize the effect of contact interactions at the unitary limit g —)■ oo, it is useful to visualize the 
coincidence manifold for three particles V^. This manifold is a structure in configuration space 
defined by the three planes qi = 92 , <?2 = 93 and qi = q^. These three planes intersect at the line 
qi — q 2 = 93 at angles of 27r/3. The point symmetry of is Hqh, which is the same as C 3 for the 
harmonic trap. See Fig. for a three-dimensional diagram of and Fig. for a two-dimensional 
diagram of relative plane cross section of V^. 

The manifold divides configuration space into six equivalent sectors, one for each order of 
particles qi < qj < qk- Denote these sectors Qs, where s is an element of S 3 expressed in permutation 
notation s = {si, 52 , 33 } and the order is 9 si < 9 s 2 < 9 s 3 - At the unitary limit, these sectors are 
dynamically isolated, but in the near-unitary limit they are connected by weak tunneling. 


4.5 Configuration Space Symmetry 

The configuration space representation 0(p) of particle permutations p S P3 map sectors Qs onto each 
other like 

Q(p)Qs = Qsp-i. (62) 

These maps are linear and continuous on Q^, mapping nearby points into nearby points in Q^, even if 
they are in different sectors. 

The set Q 3 of all possible sector exchange maps is a group isomorphic to Sg and elements are 
represented as 6 x 6 matrices M € Q 3 with a single 1 is each row and column and all the other matrix 
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Fig. 7 This figure depicts a contour plot of the relative configuration space wave function for the totally 
antisymmetric lowest energy spatial state of three particles in a harmonic well. The horizontal axis is the 
relative coordinate ri = (51 — 52)/\/2 and the vertical axis is the coordinate r 2 = (qi + (?2 — ‘2q3)/V&- The 
permutation s in curly brackets denotes each sector Qa where qa^ < qa^ < qa^- The center-of-mass degree of 
freedom is perpendicular to this plane and not depicted. The sectors are also labeled by lower-case letters for 
convenience; see Table 

elements 0. These act on the vector space of sectors like 

MQs = Qs'. 

Maps M_ that shuffle sectors are generally discontinuous transformations of Q^. The particle permuta¬ 
tions 0 {p) are the only 6 out of all 6 ! sector exchange maps that are linear and continuous on all of 
Q3. 

A particularly useful subset of sector exchange maps are ordering permutations D 3 . Instead of 
permuting particle numbers wherever they appear in the sector label s, ordering permutations permute 
the order of numbers in s, no matter what particle numbers they are. For example, the ordering 
permutation 0{AB) switches the order of the first and second number in the sector 

0(Ai?)Q{3i2} = Q{132}, 

whereas the particle permutation 0 ( 12 ) switches the positions of the numbers 1 and 2 

0(12)Q{312} = Q{321}- 

See Table for more examples of the difference between particle permutations and ordering permuta¬ 
tions. The set of ordering permutations D3 is isomorphic to S3 and acts ‘naturally’ on the sectors, i.e. 
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Table 3 This table compares the particle permutations 0{213} = 0(12) and 0(231} = 0(123) with the 
ordering permutations 0{BAC} = 0(AB) and 0{BCA} = Q{ABC). The sector notation is depicted by the 
lower-case letter in Fig.[^ Note that particle permutations have simple realizations in configuration space: two- 
cycles are reflections and three-cycles are rotations by 27r/3. Ordering permutations do not have corresponding 
orthogonal transformations of configuration space, although they form an isomorphic group to S 3 and map 
sectors onto sectors. 


Sector 

Label 

0 ( 12 ) 

0{AB) 

0(123) 

0{ABC) 

{123} 

a 

b 

b 

e 

c 

{213} 

b 

a 

a 

f 

f 

{231} 

c 

f 

d 

a 

e 

{321} 

d 

e 

c 

b 

b 

{312} 

e 

d 

f 

c 

a 

{132} 

f 

c 

e 

d 

d 


0 G D3 


0 (o)Qs = Qos- 


(63) 


Ordering transpositions {AB) and (BC) exchange sectors next to each other, i.e. the sectors most like 
to experience tunneling in the near-unitary limit. 


4.6 Kinematic Symmetry 

In the unitary limit, the six sectors Qs are dynamically isolated. There is no tunneling between them 
and is a nodal surface for every state with finite energy. More generally, because the sectors are 
dynamically isolated, the spatial Hilbert space of finite energy states /Coo C AC can be decomposed into 
sectors 


/C 


00 


0/Cs 


(64) 


where each ACg is the space of square integrable functions on Qs that vanish on two half-planes from 
V^. 
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As for two particles, a basis for each sector is provided by the snippet basis. These are restrictions 


of the totally antisymmetric states like 


to each sector 


•ims) = 


sy6<^q||^ q e Q: 

0 else 


(65) 


where tt^ is the sign of the permutation s: for two-cycles = — 1 and for three-cycles and the identity 
TTg = 1. This means that for every composition [a/SyJ of three distinct non-interacting energy eigen¬ 
states, there is a six-fold degenerate level with energy = ea + ep + e.y in the spectrum cr|°. See 

Fig. for a depiction of the relative configuration space and lowest energy level for a harmonic 
trap. 


Combining (62), (63), and (65), particle permutations p G P 3 act on the snippet basis like 


U{p) |;s 


msp 


-1 


( 66 ) 


and ordering permutations 0 € D 3 act like 


(7(0) os^ 


(67) 


From this it is clear that all p G P 3 commute with all 0 G D 3 . Therefore, every energy level 
carries a symmetry isomorphic to S 3 x S 3 = S 3 However, unlike state permutations which were a 
symmetry that depended on the compositions of the Kg irreps, the group D 3 is a symmetry of every 
finite-energy state in /Coo- 


Denote the subspace spanned by the particular snippet basis (651 by . This space carries the 


regular representation of P 3 ~ S 3 isomorphic to [3]02[21]0[1^]. Just as the non-interacting irrep space 
/C carried an additional, independent copy of S 3 called state permutation symmetry Yap-^\ that 
diagonalizes degenerate P 3 irreps, ordering permutation symmetry D 3 can be used to diagonalize the 
degenerate P 3 irreps in the snippet space . Exploiting both particle permutation symmetry and 


ordering permutation symmetry, six symmetrized snippet vectors in 


y^ G are denoted 


^5 






1 


( 68 ) 


The Young tableaux ^ filled with A, H, and C label the basis for irreps of the ordering permutation 
symmetry and the Young tableaux Y filled with particle numbers label the basis for irreps of the particle 
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Fig. 8 This figure depicts a contour plot of the relative configuration space wave function for the lowest 
degenerate states in the unitary limit of the contact interaction for three particles in a harmonic well. The 
horizontal axis is the relative coordinate ri and the vertical axis is the coordinate r 2 . The states are depicted 
in the same order as ( |68[ ). 

permutation symmetry. The exact expression for these symmetrized vectors in terms of snippet bases 
vectors depends on the choice of subgroup chains for the complete set of commuting observables. For 
symmetric wells a good choice of observables is , and = U(AC) because within a 

snippet subspace with composition [a/JyJ total parity inversion U is proportional to U{AC) 


fl 112) r) = -TT^irpn^UiAC) ||; 2) r) . 


(69) 


With this convention, the states 




, and 


have the opposite parity from the 


totally antisymmetric state 




Hy See Fig. 8 for a depiction of the six states in using the 

complete set of commuting observables to diagonalize the lowest energy level 

in the harmonic well. 


Note that one-component bosons can only populate the state 


M \a\b\c\ \ and one-component fermions 


mh). To include spin, consider spin-1/2 distinguishable particles with 


can only populate the state 
5 = C®. Each energy level X 5 is 48-fold degenerate because there are eight spin states for every 


spatial state. Following the method described for the non-interacting, identical particles with spin in 
Section 4.2, using the Clebsch-Gordan series for S 3 one can show that there are four two-component 
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bosonic (fermionic) states with spatial symmetry [3] ([1^]) and total spin s = 3/2 and four with spatial 
symmetry [21] and total spin s = 1/2. 


4.7 Near-Unitary Limit 


Another use for ordering permutation symmetry is looking at the near-unitary limit. In this limit 
there is tunneling between adjacent sectors. For an asymmetric well, the tunneling operator T can be 
parameterized 

f=-tU{AB)-uU{BC)-{t + u)U{e), (70) 


where t and u are the tunneling amplitudes for the first and second particle and the second and third 
particle [29]: 



9 


6 

u = - 
9 








2 


91 =® 

2 


92=93 


(71) 


The coefficients t and u depend on trap shape through the one-particle eigenstate wave functions, and 


for a symmetric trap t = u. The third term in (70) which is proportional to the identity renormal¬ 
izes the energy shift such that the totally antisymmetric state undergoes no change in energy. The 
first order energy shifts are given by the eigenvalues of this matrix formed by the matrix elements 


l;2}F 


T 


5 F). The eigenvalues and degeneracies are 


Z\A[ 3 ] = —2t — 2u (singly-degenerate), 

4\A[2i] = —t — u — \/t'^ — tu + v? (doubly-degenerate), 

AE[ 21 ]' = —t — u + \/t'^ — tu + v? (doubly-degenerate), 

Z\A[i 3] = 0 (singly-degenerate). (72) 


The first and last eigenvalues correspond to symmetric and antisymmetric states, and the middle 
two eigenvalues are for two different two-dimensional mixed-symmetry irreps. Since the ratio of t to 
u depends on trap shape, the first-order energy depends on trap shape for asymmetric traps. The 
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symmetric [3] and antisymmetric states [1^] are universal for any trap, but energy eigenstates with 
mixed symmetry [21] are not, although there are universal algebraic expressions in terms of t and u. 
However, for symmetric traps t = u and the eigenvalues become 

Z\if[ 3 ]± = —At (singly-degenerate), 

= —3t (doubly-degenerate), 

[ 21 ]± = —t (doubly-degenerate), 

Z\i?[i3]T = 0 (singly-degenerate). (73) 

With the addition of parity symmetry, we find that although the eigenvalues depend on the trap 
shape through t, the states do not. The ordering permutation observable U{AC) ^ 77 is sufficient to 
distinguish the two copies of the [21] irreps. The energy shifts —3t and 0 are for states with the same 

and the energy shifts —At and —t are states with 
opposite parity -■Ka'Kp'K^. These results agree with recent work on the near unitary limit for harmonic 
traps [ini 1301 EZl EH] where t = u = 3^/(2^v^ff)- 

5 Conclusion of Part I 

What lessons have been learned from this analysis of the symmetries for one, two and three trapped 
particles? What methods can be extended to more particles, and what methods lose utility? 

One set of observations is about how symmetries of the individual particles build into symmetries of 
the total system. Any system with a Hamiltonian that is a sum of identical sub-Hamiltonians, each with 
symmetry Gi, must have at least the minimal the symmetry structure Pjv x (Gi)^^. In particular, the 
configuration space symmetry and the kinematic symmetries of two or three non-interacting particles 
satisfy 

3 Pjv X Cf ^ and 3 Pw x Kf (74) 

When these groups are greater than the minimal symmetry, then that signals the emergence of a true 
multiparticle symmetry. For two and three particle systems, the groups 0° and Cg are finite order 
point groups in two and three dimensions, realized by orthogonal transformations of configuration 
space. Unless there are emergent symmetries or accidental degeneracies, the irreps of K 2 and Kg 


parity TTaTTpTT~^ as the original fermionic state 
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are sufficient to explain the degeneracy pattern of the non-interacting energies. The main difference 
between the case of two and three particles is that for three particles there can be multiple irreps of 
the symmetric group S 3 with the same energy, even without emergent or accidental symmetries. State 
permutation symmetry is introduced to build observables that distinguish these degenerate spaces. 
All of this is extended to more particles in the sequel article using the mathematics of permutation 


modules, and the minimal construction (74) is shown to be algebraically universal, i.e. a complete set 


of observables that is independent of trap shape can always be built out of single particle operators. 
For systems with Galilean-invariant, spin-independent two-body interactions, then the symmetry 


(74) is broken and algebraic universality is lost. However, the symmetry that remains is always at least 


Cat 3 Pat x Ci and Ka- D Pjy x Ki. 


(75) 


The configuration space symmetries again have geometrical realizations as point groups that preserve 
both the trap equipotentials and the coincidence manifold. The irreps of the kinematic symmetries are 
irreps of Sa with additional quantum numbers inherited from the single-particle trap, e.g. parity for 
symmetric traps and center-of-mass quantum number for harmonic traps. The energy eigenvalues and 
eigenstates are certainly not universal, but the symmetries can be used to decrease the computational 
scale of exact diagonalization or perturbation theory by using basis vectors from the non-interacting 
spatial Hilbert space reduced into irreps of Ka . 

For weak interactions, two and three particles have different degrees of algebraic universality. The 
two-body matrix elements, which depend on the trap and the interaction, are required in order to find 
the specific energy shift at first order. However, the splitting pattern and the energy eigenstates are 
independent of interaction and trap shape for two particles and are therefore algebraically universal 
in a strong sense because the two-body matrix elements are not required. For three particles, there 
is less universality under weak perturbations. Some splitting features remain independent of trap 
and interaction, but the full solution for some first-order eigenstates requires algebraic expressions 
involving the two-body matrix elements. In the sequel article even this very limited nniversal feature 
is hypothesized to fail for > 5. 

The contact interaction restores some algebraic universality, and this can be partially understood 
as a manifestation of the state permutation symmetry of the two-body matrix elements. Also, for the 
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contact interaction the unitary limit is sensible because the delta interaction enforces nodal surfaces in 
configuration space. These surfaces prevent particles tunneling past each other, and break configuration 
space into sectors with a particular order. The emergent combination of ordering permutation symmetry 
and particle permutation symmetry provides enough structure for complete algebraic solution of the 
Hamiltonian in the algebraic limit assuming knowledge of the single particle spectrum. The near unitary 
mapping is also algebraically universal for two particles, but for three particles in an asymmetric trap 
information about the trap shape (or equivalently the single particle wave functions) is required to 
solve for the level splitting. For symmetric traps, the extra information provided by parity is enough 
to restore algebraic universality. 
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